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1. Introduction 



An important open problem of the AdS / CFT correspondence is to understand the 
finite-size spectrum of the AdS 5 x S 5 superstring. Recently, there has been further 
significant progress in this direction. First, the four-loop anomalous dimension of 
the Konishi operator was computed by means of generalized Liischer's formulae 
[0, f|, (see also ||-[[nj for other applications of Liischer's approach), and the result 
exhibits a stunning agreement with a direct field-theoretic computation |L6], |T7|j . 
Second, the groundwork for constructing the Thermodynamic Bethe Ansatz (TBA) 
|P-8|| , which encodes the finite-size spectrum for all values of the 't Hooft coupling, 
has been laid down, based on the mirror theory approach 1 ||20|| . Most importantly, 
the string hypothesis for the mirror model was formulated pll] an d used to derive 



TBA equations for the ground state [22]- [25]. Also, a corresponding Y-system |p6| 
was conjectured p7j, and its general solution was obtained [p8| |. The AdS/CFT 
Y-system has unusual properties, and, in particular, is defined on an infinite-genus 
Riemann surface [^, . 

The derivation of the TBA equations is not yet complete, because the equations 
pertain only to the ground state energy (or Witten's index in the case of periodic 
fermions) and do not capture the energies of excited states. Therefore, one has to find 
a generalization of the TBA equations that can account for the complete spectrum 
of the string sigma model, including all excited states. 

Here we continue to explore the mirror TBA approach. In particular, we will 
be interested in finding the TBA integral equations which describe the spectrum 
of string states in the sl(2) sector. An attempt in this direction has been already 
undertaken in [25] and the emerging integral equations have been used for numerical 
computation of the anomalous dimension of the Konishi operator |30 |. However, the 
subleading term in the strong coupling expansion in this result disagrees with the 
result by [[3TJ obtained by string theory means. There exists yet another prediction 
for this subleading term, which differs from both [RO] and pi|. All these results 



are based on certain assumptions which require further justification. This makes 
urgent to carefully analyze the issue of the TBA equations for excited states, and to 
better understand what happens on the string theory side. 

In this paper we analyze two-particle states in the sl(2) sector. First, we show 
that each state is governed by its own set of the TBA equations. Second, we provide 
evidence that for each state there are infinitely-many critical values of 't Hooft cou- 
pling constant A, and that the excited states integral equations have to be modified 
each time one of these critical values is crossed. 2 Performing careful analysis of two- 



1 The TBA approach in the AdS/CFT spectral problem was advocated in [|19| where it was used 
to explain wrapping effects in gauge theory. 

2 Existence of such critical values was observed in the excited-state TBA equations for perturbed 
minimal models p%. We thank Patrick Dorey for this comment. 
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particle states in a region between any two neighboring critical points, we propose 
the corresponding integral equations. 

The problem of finite-size spectrum of two-dimensional integrable models has 



been studied in many works, see e.g. 33 - 49 . To explain our findings, we start 



with recalling that for some integrable models the inclusion of excited states in the 
framework of the TBA approach has been achieved by applying a certain analytic 
continuation procedure [37| . This can be understood from the fact that the 
convolution terms entering the integral TBA equations exhibit a singular behavior 
in the complex rapidity plane, the structure of these singularities does depend on the 
value of the coupling constant. This leads to a modification of the ground-state TBA 
equations, which, indeed, describe the profile and energies of excited states. Here we 
intend a similar strategy for the string sigma model. 

To derive the TBA equations for excited states, we propose to use a contour 
deformation trick. In other words, we assume that the TBA equations for excited 
states have the same form as those for the ground state with the only exception 
that the integration contour in the convolution terms is different. Returning the 
contour back to the real rapidity line of the mirror theory, one picks up singularities 
of the convolution terms which leads to modification of the final equations. The 
original contour should be drawn in such a way, that the arising TBA equations 
would reproduce the large L asymptotic solution (where L is the size of the system). 

Recall that the TBA equations for the string mirror model |22| are written in 
terms of the following Y-functions: Ig-functions associated with Q-particle bound 
states, auxiliary functions Yq\ vw for Q|uw-strings, Yq\ w for Q|u>-strings, and Y± for 
2/±-particles. The Y-functions depend on the 't Hooft coupling A related to the string 
tension g as A = 4n 2 g 2 . As we will see, the analytic structure of these Y-functions 
depends on g and plays a crucial role in obtaining the TBA equations for excited 
states. 

Most conveniently, the large L asymptotic solution for the Y-functions is written 
in terms of certain transfer-matrices associated with an underlying symmetry group 
of the model In the context of the string sigma model the corresponding 

asymptotic solution was presented in P7| . We will use this solution to check the 
validity of our TBA equations. 

Our analysis starts from describing physical two-particle states in the $1(2) sector. 
It appears that for the $1(2) states the functions Yq\ vw (u) play the primary role in 
formulating the excited state TBA equations. Analyzing the asymptotic solution, 
we find that each Yq\ vw has four zeroes in the complex w-plane. With g changing, 
the zeroes change their position as well, and at certain critical values g = g cr they 
give rise to new singularities in the TBA equations which resolution results in the 
appearance of new driving terms. The critical values g cr are defined as values of g 
at which Yq\ vw (u) acquires two zeros at u — ±i/g cr : Yqi vw (±i/ g cr ) = 0. 
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We show that at weak coupling g ~ all two-particle states can be organized 
into an infinite tower of classes k — I , II , . . . ,00, see Table 1. 



Type of a state 


Y-functions 


Number of zeros 


I 




2 


II 




2+2 


III 


^l|mLi> ^2\vw7 ^3 \vw 


4+2+2 


IV 


Y Y Y V 
1 l\vwi ± 2\vun J 3t>uij 1 4\vw 


4+4+2+2 








k — > 00 


Y V i 

1 l\vwi 1 2\vwi • • • 


4+4+ . . . 



Table 1. Classification of two-particle states in the s[(2)-sector at g ~ 0. The 
right column shows the number of zeros which the corresponding asymptotic 
Y^l^-functions from the middle column have in the physical strip |Im(u)| < 
1/g. States of type I are called "Konishi-like" . States of type II, III, IV, . . . 
correspond to larger value of k, see section 3. 



Each class is unambiguously determined by number of zeroes of YM|?™-functions in 
the strip |Im(w)| < 1/g. In particular, for states of type I only Fi^-function has two 
zeroes in the physical strip. We call all these states "Konishi-like" because they share 
this property with the particular string state corresponding to the Konishi operator. 

Our results disagree with that by p5], [3(J in the following two aspects: First, 
the integral equations for excited states in the 5 1(2) sector do not have a universal 
form, even for two-particle states. Second, we face the issue of critical points. When 
a critical point is crossed, the compatibility of the asymptotic solution with the 
integral TBA equations requires modification of the latter. The equations proposed 
in 



25| capture only type I states and only below the first critical point. 



To find approximate locations of the critical values, we first solve the asymptotic 
Bethe-Yang equations |50| (which include the BES/BHL dressing phase |51~| , |52j| ) for 



some states numerically from weak to strong coupling and obtain the corresponding 
interpolating curve u = ui(g), where u\ is the rapidity of an excited particle in 
string theory. The second particle has rapidity 112 = —U\ due to the level matching 
condition. Further, we compute the Y-functions on the large L asymptotic solution 
corresponding to this two-particle excited state and study their analytic properties 
considered as functions of g, and, in particular, determine approximately the critical 
values. 

We also note that with the coupling increasing more and more critical points 
get crossed which leads to accumulation of zeroes of Fm|™'s in the physical strip. 
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Apparently, as the asymptotic solution indicates, when g tends to infinity the zeroes 
move towards the points ±2, so that the latter points behave as an attractor for 
zeroes of all l^i^-functions. 

Estimation based on the large L asymptotic solution gives A ~ 774 for the first 
critical value corresponding to the Konishi operator. In the weak-coupling region 
below the first critical point, the integral equations for Konishi operator we obtain 
seem to agree with that of |25]. However, these weak-coupling equations become 
inconsistent with the known large L asymptotic solution once the first critical point 
is crossed, and have to be modified. Consequently, the derivation of the anomalous 
dimension for the Konishi operator at strong coupling requires re-examination. Of 
course, the existence of critical points is not expected to violate analyticity of the 
energy E(X) of a string state considered as the function of A, but it poses a question 
about the precise analytic behavior of -E'(A) in the vicinity of critical points. 

We discuss both the canonical and simplified TBA equations. The canonical 
equations P2L E3, E31 follow from the string hypothesis for the mirror model [ETJ by 



using the standard procedure, see e.g. H53| . The simplified equations [|22], ^3) obtained 
from the canonical ones have more close relation to the Y-system. It turns out that 
the simplified equations are sensitive only to the critical points defined above. In 
contrast, the canonical equations have to be modified when crossing not only a critical 
point but also what we call a subcritical point g cr . A subcritical point g cr is defined 
as the value of g at which the function Yq\ vw (u) acquires zero at u = 0. Hence, 
in comparison to the canonical equations, the simplified equations exhibit a more 
transparent analytic structure. In addition to locality, this is yet another reason why 
we attribute to the simplified equations a primary importance and carry out their 
analysis in the main text. To study the exact Bethe equations which determine the 
exact, i.e. non-asymptotic, location of the Bethe roots, we find it advantageous to 
use a so-called hybrid form of the TBA equations for Yg-functions. This form is 
obtained by exploiting both the canonical and simplified TBA equations. 

Recently, the finite-gap solutions of semi-classical string theory have been nicely 
derived |54[] from the TBA equations [25]. This raises a question why modifications 
of the TBA equations we find in this paper were not relevant for this derivation. 
We have not studied this question thoroughly. However, one can immediately see 
that there is a principle difference between states with finite number of particles and 
semiclassical states composed of infinitely many particles. Namely, at strong coupling 
the rapidities of two-particle states fall inside the interval [—2, 2], while those of semi- 
classical states are always outside this interval. Thus, the modification of the TBA 
equations discussed in this paper might not be necessary for semi-classical states. It 
would be important to better understand this issue. 

The paper is organized as follows. In the next section we explain our criteria for 
a choice of the integration contour in the excited states TBA equations. In section 3 
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we discuss two-particle states in the $1(2) sector and the corresponding asymptotic 
Y-functions. By analysing analytic properties of the Y-functions, we determine the 
critical and subcritical values of the coupling constant both for the Konishi and 
for some other states. In section 4 we present the simplified TBA equations for 
Konishi-like states and we explain why and how their form depends on the value of 
the coupling constant. In section 5 we generalize this discussion to arbitrary two- 
particle states from the s 1(2) sector. In section 6 we summarize the most essential 
properties of the AdS/CFT Y-system implied by the TBA equations under study. 
Finally, in Conclusions we mention some interesting open problems. The definitions, 
treatment of canonical equations, and further technical details are relegated to eight 
appendices. 



2. Contour deformation trick 

The TBA equations for the AdSs x S 5 mirror model are written for Y-functions which 
depend on the real momentum of the mirror model. The energy of string excited 
states obviously depends on real momenta of string theory particles, and to formu- 
late the TBA equations for excited states one also needs to continue analytically the 
Y-functions to the string theory kinematic region. To visualize the analytic continu- 
ation it is convenient to use the ZQ-tori because the kinematic regions of the mirror 
and string theory Q-particle bound states (Q-particles for short) are subregions of 
the z-torus, see Figure 1. In addition, the Q-particle energy, and many of the kernels 
appearing in the set of TBA equations are meromorphic functions on the correspond- 
ing torus. The mirror Q-particle region can be mapped onto a u-plane with the cuts 
running from the points ±2 ± -Q to ±oo, and the string Q-particle region can be 
mapped onto a w*-plane with the cuts connecting the points —2 ± and 2 ± ^Q, 
see Figure 1. Since the cut structure on the planes is different for each Y-function, 
they cannot be considered as different sheets of a Riemann surface. The z^-torus 
can be glued either from four mirror w-planes or four string w*-planes. 

As was shown in [pOj] , z-torus variables corresponding to real momenta of the 
mirror and string theory are related to each other by the shift by a quarter of the 
imaginary period of the torus 

z = z* + y, (2.1) 

where z is the variable parametrizing the real momenta of the mirror theory, and z* is 
the variable parametrizing the real momenta of the string theory. The line (— oo, oo) 
in the string u*-plane is mapped to the interval Re(z*) G (— ^-), Im(z*) =const 
on the z-torus, and we choose z* in the string region to be real. Then, the interval 
Re(z) G {—^2i ^r), Im(,2) = || of the mirror region is mapped onto the real line of 
the mirror w-plane. 
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Figure 1: These are the mirror and string regions on the z-torus. They are in one-to-one 
correspondence with the it-planes. The boundaries of the regions are mapped to the cuts. 



It is argued in |36|, [3^] that the TBA equations for excited states can be obtained 



from the ones for the ground state by analytically continuing in the coupling constants 
and picking up the singularity of proper convolution terms. We prefer however to 
employ a slightly different procedure which we refer to as the contour deformation 



trick. We believe it is equivalent to [Bq, 37 1. It is based on the following assumptions 



The form of TBA equations for any excited state and the expression for the 
energy are universal. TBA equations for excited states differ from each other 
only by a choice of integration contours of convolution terms and the length 
parameter L which depend on a state. 



• The choice of the integration contours and L is fixed by requiring that the large 
L solution of the excited state TBA equations be given by the generalized 
Luscher formulae, that is all the Y-functions can be written in terms of the 
eigenvalues of the transfer matrices. The integration contour depends on the 
excited state under consideration, and in general on the values of 't Hooft's 
coupling and L. 

• An excited state is completely characterized by the five charges it carries and a 
set of N real numbers pk which are in one-to-one correspondence with momenta 
p° k of iV Q-particles in the small coupling limit g — > 0. The momenta p° k are 
found by using the one-loop Bethe equations for fundamental particles and 
their bound states. For finite values of g the set of pu is determined by exact 
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Bethe equations which state that at any p = pk the corresponding ^-functions 
are equal to —1. 

We consider only the simplest case of two-particle excited states in the sl(2) 
sector of the string theory because there are no bound states in this sector and the 
complete two-particle spectrum can be readily classified. The physical states satisfy 
the level-matching condition which for two-particle states takes a very simple form: 
Pi = —P2, or u\ = —U2 or z*i = —z* 2 , depending on the coordinates employed. The 
TBA equations we propose in next sections are valid only for physical states. 



3. States and Y-functions in the sl(2) sector 

To fix the integration contour one should choose a state and analyze the analytic 
structure of the large L Y-functions which we refer to as the asymptotic Y-functions. 
We begin with a short discussion of two-particle states in the sl(2) sector. 

3.1 Bethe- Yang equations for the sl(2) sector 

There is only a single Bethe- Yang (BY) equation in the s[(2)-sector for two-particle 
physical configurations satisfying the vanishing total momentum condition p±+p 2 = 
that can be written in the form [[55 



P J S. m {p,-p) e^ J+1 > 



1+ 1 



.+ 2 



1 + 



'—<r(p, -pY 



(3.1) 



where p = p\ > 0, J is the charge carried by the state, a is the dressing factor, and 
xf are defined in appendix IO. Taking the logarithm of the equation, one gets 



(3.2) 



1 + ^ 

ip(J + 1) — log -j 2% 9(p, —p) = 2iri n . 

1 + — 0. 



where 6 = - log a is the dressing phase, and n is a positive integer because we have 
assumed p to be positive. As was shown in [ 56) , at large values of g the integer n is 
equal to the string level of the state. 

As is well known, in the small g limit the equation has the obvious solution 



Pi 



2nn 

7+i 



n 



J + l 



(3.3) 



where [x] denotes the integer part of x, and the range of n is bounded because the 
momentum p can only take values from to tt. The corresponding rapidity variable 
uj )Tl in the small g limit takes the following form 



u J,n ~^ ~^ u J,n 



u 



J.n 



cot ■ 



nn 



J+l 



(3.4) 
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Figure 2: These are the plots of u which solves the BY equation for the Konishi state. 



Thus, any two-particle state in the sl(2) sector is completely characterized by the 
two integers J and n. In particular, in the simplest J = 2 case corresponding to a 
descendent of the Konishi state n can take only one value n = 1, and the small g 
solution is 

p£,i = y> ^> i = 7r (3 - 5) 

The BY equation (|3.2| ) can be easily solved perturbatively up to any desired order 
in g, and numerically up to very large values of g. We have used the BES series 



representation [51]] for the dressing phase for perturbative computations, and the 
DHM integral representation |J7| for the numerical ones 3 . For the Konishi state, the 
perturbative solution up to g 16 -th can be found in appendix K2. 



We have solved numerically the BY equation for the Konishi state for i < g < 
1000 with the step i for i < g < 10, the step 1 for 10 < g < 100, and the step 
10 for 100 < g < 1000. In Figure 2 we show the results up to g = 100. For greater 
values of g nothing interesting happens, and the solution can be approximated by 
asymptotic formulae from |56|, [32] |5J|, see appendix |8.2| for more details. We then 
applied the Interpolation function in Mathematica to have U2 t i as a smooth function 
of g. Using the function, one can find that u 2 ,i(g) decreases up to g ~ 0.971623, and 
then begins to increase and at large g it asymptotes to u — 2. 

The functions uj )U (g) for other values of J and n have similar g dependence. The 
only exception is the case J = 2n — 1 where the exact solution of the BY equation 
is P2n-i,n = 7T and, therefore, U2 n -i,n = 0. The TBA equations we propose below are 
not in fact valid for these (2n — 1, n) states. 

The perturbative and numerical solutions for u j, n (g) can be easily used to analyze 
the behavior of Y-functions considered as functions of g. In particular it is easy to 
determine if some of them become negative for large enough values of g. 



^Thc DHM representation can be also readily used for perturbative computations. 
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3.2 Y-functions in the sl(2) sector 

Let us recall that the TBA equations for the A0IS5 x S 5 mirror model involve Yq- 
functions for momentum carrying Q-particle bound states, and auxiliary functions 
Yq\vw ^ or Ql^w-strings, Yq?^ for Q|w;-strmgs, and for y±-particles. The index 
a = 1, 2 reflects two su(2|2) algebras in the symmetry algebra of the light-cone string 



sigma model. The TBA equations [52j depend also on the parameters h a which take 
care of the periodicity condition of the fermions of the model . For the s[(2) sector 
the fermions are periodic, and from the very beginning one can set the parameters h a 
to because there is no singularity at h a = in the excited states TBA equations. 

For the s[(2) states there is a symmetry between the left and right su(2|2) aux- 
iliary roots, and, therefore, all Y-functions satisfy the condition 

Y v (1) = Y v (2) = Y v , (3.6) 

where V denotes a Y-function of any kind. 

The string theory spectrum in the s[(2) sector is characterized by a set of N real 
numbers z*k or u k corresponding to momenta of N fundamental particles in the limit 
g — > 0. According to the discussion above, these numbers are determined from the 
exact Bethe equations 

Yi(z* fe ) = Y u (u k ) = -1, k = l,...,N, (3.7) 

where Yi(z) is the Y-function of fundamental mirror particles considered as a function 
on the z-torus, and Yi„ (u) denotes the Yi-function analytically continued to the string 
w-plane. 

In the large J limit the exact Bethe equations must reduce to the BY equations, 
and it is indeed so because the asymptotic s((2) Yg-functions can be written in terms 
of the transfer matrices defined in appendix |S.3| as follows [[| 

rp / I ->\ 2 _ N 

W = e- J ^ W n ;% ( y : = e-^T Q Muf \{S^u t ) , (3.8) 
n,;=i^(2)(^^) f=i 

where v is the rapidity variable of the mirror w-plane and Sq is the energy of a 
mirror Q-particle. S^S, denotes the S-matrix with the first and second arguments 
in the string and mirror regions, respectively. Tq^{v\u) is up to a factor the trace of 
the S-matrix describing the scattering on these string theory particles with a mirror 
Q-particle or in other words the eigenvalue of the corresponding transfer matrix. 
The BY equations then follow from the fact that £u(v,k) = —ipk, and the following 
normalization of T^i 

N 

T 1A (u tk \u) = l =^ -l = e iJ »'Y[Sl$(u k ,u i ), (3.9) 
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where = Uk, and the star just indicates that one analytically continues T^i to the 
string region. Then, S s Z 2 ?(uk,Ui) = S S [(2)(uk,Ui) is the usual sl(2) sector S-matrix 
used in the previous subsection. 

Let us also mention that Tq i has the following large v asymptotics 

T q>1 {v\u)^q(i-Y[^\ , v^oo, (3.10) 

and therefore it goes to if the level-matching is satisfied. 

Then, in the large J limit all auxiliary asymptotic sl(2) Y-functions can be 



written in terms of the transfer matrices as follows 127 



Ti rri rri rri rri rri rri 

yo = 2^ y° = 3 Y° -= Q+2 ' 1 QA Y° — 1 ' Q+21 ' Q 



Ti J ~ + rp rp 1 Q\vW rp 7 - IJ\W rp rp 

1,2 i 3,2 J l,2 " J-Q+1,2 ' I 2,Q+1 I 0,Q+1 

The transfer matrices T a ^ s can be computed in terms of T flj i by using the Bazhanov- 
Reshetikhin formula |JT|, see appendix ^]3]for all the necessary explicit formulae. 

An important property of the Y-functions for vw- and ^-strings is that they 
approach their vacuum values as v — > oo 

M M 

Y M \vw(v) ^M(M + 2), Y M Uv) ->■ M(M + 2) , v -> oo , < Jm(v) < — . 

9 9 

Now we are ready to analyze the dependence of asymptotic Y-functions on g. 
Recall that they depend on the rapidities Uk which are solutions of the BY equations. 

We begin with the small g limit where the effective length goes to infinity, and 
one can in fact trust all the asymptotic formulae. It is convenient to rescale v and Uk 
variables as v — >■ v/g, Uk — > Uk/9 because the rescaled variables are finite in this limit. 
Let k = U\ — —U2 be the rescaled rapidity of a fundamental particle. According to 
the previous subsection, in the small g limit they are given by Uj n , eq.(|3.4j). 

The most important functions in the $1(2) case are Ym\ V w, and we find that for 
N = 2 they exhibit the following small g behavior in the strip — M < lm(v) < M 



[M 2 — l + v 2 — k 2 ] [(M + 2) 2 - 1 + v 



2 -K 2 



Y M{vM = M(M + 2) L_ _ _ _ _ _ ^ _ j_ _ _ _ _ _i + W . (3 . n) 

The leading term has the correct large u-asymptotics and four apparent zeros at 

v = ±Vk 2 - M 2 + 1 , v = ±a/ k 2 — (M + 2) 2 + 1 . 

One can see that Yi| utu -function always has at least two real zeros at v = ±k. Other 
zeros of Ym V w-t unctions can be either real or purely imaginary depending on the 
values of M and k. It appears that the form of simplified TBA equations depends 
on the imaginary part of these zeros, and we will see in next sections that if a pair 
of zeros v = ±r fall in the strip |Im(r)| < 1 then the equations should be modified. 
Thus, we are lead to consider the following three possibilities 
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If M 2 — 2<k 2 < (M + 2) 2 -2 then Y M{vw has two zeros at v = ±aA 2 - M 2 + 1 
that are in the strip \lm(v ) | < 1. In terms of the integers J and m characterizing 
two-particle states one gets the condition 



VM 2 -2 < cot < v / (M + 2) 2 -2 . (3.12) 



2. If k 2 < M 2 — 2 -<=>- cot j^j < \/M 2 — 2 then Y M \ VW does not have any zeros 
in the strip |Im(t>)| < 1. 

3. If k 2 > (M + 2) 2 — 2 cot > v/(M + 2) 2 - 2 then y M |™ has four zeros 
in the strip |Im(i>)| < 1. 

Some of these zeros can be real, and in fact the canonical TBA equations take 
different forms depending on whether the roots are real or imaginary. 

Classification of two-particle states at g ~ is presented in Table 1. The type of 
a state is determined by how many zeroes of >M|i™-functions occur in the physical 
strip and it depends on J and n. 

Consider a two-particle state with k = u° Jn for some (J,n). Table 1 shows that 
there exists a number m > 1, equal to the maximal value of M the condition ( [3.12] ) 
is satisfied. Then both Y m \ vw and Y m -\\ vw have two zeros, all Y k \ vw with k < m — 2 
have four zeros, and all Yu^-functions with k > m + 1 have no zeros in the strip 
|Im(u)| < 1. For example, among the states with (J,n = 1) at small coupling, the 
states of type I are found if and only if J < 4. The type II is found for 5 < J < 7, 
type III for 8 < J < 11, and type IV for 12 < J < 14. In particular, Y^ vw for the 
state (8, 1) has two real zeros and two imaginary zeros in the strip |Im(t>)| < 1, and 
Y\\ vw for the state (J > 9, 1) has four real zeros. 

As for the Konishi state with J = 2 and n — 1 only Yi^-function has two 
zeros and all the other y^^to-functions have no zeros at small coupling. Let us also 
mention that at g — the l^to-function of the state (5, 1) (and in general of any 
state (6k — 1, k)) has a double zero at v = 0. This double zero however is an artifact 
of the perturbative expansion, and in reality Y 2 \ vw has two imaginary zeros for small 
values of g equal to ~ ±ig\/3. For the state with J = 6 and n = 1 both Yx\ vw and 
Y 2 \ vw have two real zeros. 

3.3 Critical values of g 
Evolution of zeros 

Now we would like to understand what happens with Y M \ vw -iwactioias when one starts 
increasing g. To this end one should use numerical solutions of the BY equations 
discussed at the beginning of this section. We also switch back to the original u 
variables because they are more convenient for general values of g, and refer to the 
strip |Im(u)| < 1/g as the physical one. 
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We find that for finite g any Y^i ^-function has four zeros which are either real 
or purely imaginary. We could not find any other complex zeros. The four zeros of 
Y k \ vw are split into two pairs, and the two zeros in a pair have opposite signs, and are 
either real or complex conjugate to each other. We denote the four zeroes of Y k \ vw 

by r^p and ff\ where j — 1,2 and k — 1 , 2 , If the four zeros are real then 

are the zeros of Y k \ vw which have a larger absolute value than fj k \ If only two zeros 
are real then we denote them as and the imaginary zeros as fj k \ Finally, if the 
four zeros are imaginary then are the ones closer to the real line than the second 
pair Vj . 

The locations of the zeros depend on g, and we should distinguish two cases. 
We observe first that if two zeros are real at g ~ then they are of order 1/g, and 
obviously outside the interval [—2,2]. With g increasing they starts moving toward 
the origin, and at some value of g they reach their closest position to the origin which 
is inside the interval [—2,2]. Then, for larger g they remain inside the interval but 
begin to move to its boundaries and reach them at g — oo. In the second case, one 
considers a pair of imaginary zeros at g ~ 0. With g increasing they start moving 
toward the real line, and at some value of g they get to the origin and become a 
double zero. Then, for larger g they split and begin to move to the boundaries of the 
interval [—2, 2], and reach them at g — oo. The only exception from this behavior we 
find is the (^-dependence of the two zeros of Y"2|m«-function for the states (6k — 1, k) 
that are equal to ±i\/3 at small g. These zeros become real at very small value of g. 
Then they start moving to ±2, cross the boundaries of the interval [—2, 2], and reach 
their maximum. After that they behave as zeros of all the other Y^^-functions. 
Thus, at very large values of g all zeros of any Y k \ vw - function are real, inside the 
interval [—2,2], and very close to ±2. 

The pairs of the zeros of different Y^my-functions are not independent, and 
satisfy the following relations 

rf-V = r f +1) , fc = 2,...,oo. (3.13) 

Therefore, the zeros of Y M \ vw -i unctions can be written as follows 

{r (D , r (3) } . {r (2) > r (4) } . . {r (*-l) r ( fe+ D } . {r « r ( fc+2)} . {r ( fc +D r (*+3) } . (314) 

so that Y k \ vw has the zeros {r^p ,rj k+2 ^}. 

These zeros have a natural ordering. If we assume for definiteness that the zeros 
with j = 1 have negative real or imaginary parts, then they are ordered as 

r W -< r f0 -< r W -< . . . -< r W -< r f +1) -<-.., (3.15) 

where r[ k) -< rf +1) if either Re(rf ) ) < Re(rf +1) ) or Im(rf°) > Im(rf +1) ). It is 
important that the zeros never change the ordering they have at g ~ 0. In particular, 
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Yi\ vw always has two real zeros rj = Uj which are Bethe roots. They are the largest 
(in magnitude) zeros among all Y^^-functions, and are the closest ones to ±2 at 
large g. 



Initial condition — > 


Y V 

1 l\vw i 1 2\vw 


2+2 




y y y 
1 X\vwi 1 2\vwi 1 3\vw 


4+2+2 


9 I 


Y Y Y Y i 

1 \\vw i 1 2 \vw ) ± 3\vwi 1 4\vw 


4+4+2+2 




Y\ \vwj ^2\vwi ^3|t)ui> ^4|mo> ^5™ 


4+4+4+2+2 




Y Y 

1 \\vw ) 1 2\vw ) • ■ ■ 


4+4+ . . . 



Table 2. Evolution of a two-particle states in the sl(2)-sector with respect 
to g. At g ~ a state has a certain number of Yj^| ^-functions with zeroes 
in the physical strip. Increasing the coupling, the critical points get crossed 
which leads to accumulation of zeroes of Ym\vw ,s i n the physical strip. This 
phenomenon can be called "Y-function democracy" . 

In addition we find that the functions below have either zeros or equal to —1 at 



locations related to r 



(A:) 



Y, 



k\vw 



(fc+1) 



9 



n +1 (rf +1) ) 



Y±(r^ 



3 > 



0.(3.16) 



As will be discussed in the next section the equations Yk\ vw {r^ ± ^) = — 1 lead to 
integral equations which play the same role as the exact Bethe equations Yi(%) = — 1 
and allow one to find the exact location of the roots r^ 1 '. 

Let us finally mention that nothing special happens to YM|™-functions. 

Critical and subcritical values 



Let m again be the maximum value of M the condition ( |3.12| ) is satisfied. According 
to the discussion above for any two-particle state there is a critical value of g such 
that the function Y m+ i\ vw which had no zeros in the physical strip for small values 
of g, acquires two zeros at u = ±i/g. At the same time Y m _^ vw also acquires zeros 
at u = ±i/g. At a slightly larger value of g the two zeros that were at u = ±i/g 
collide at the origin, and Y m+ i\ vw and 5^-11^ acquire double zeros at u — 0. Then, 
the double zeros split, and both Y m \ vw and Y m+ i\ vw have two real zeros, and Y m _i\ vw 
has four. Increasing g more, one reaches the second critical value of g such that the 
functions 3^+21^ an d Y m \ vw acquire zeros at u — ±i/g, see Table 2. 

This pattern repeats itself, and there are infinitely many critical values of g which 
we denote as g r j™ and define as the smallest value of g such that for a symmetric 
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configuration of Bethe roots the function Y m+r \ vw acquires two zeros at u — ±i/g. 
The subscript J,n denotes a state in the s((2) sector, and they determine m. 

The critical values of g can be also determined from the requirement that at g = 
g r j™ the function l + Y m+r _^ vw has a double zero at u = 0: 1 +Y m+r ^x\ vw (0, g r j™) = 0. 
This condition is particularly useful because the value of the Y-functions at u = 
can be found from the TBA equations, see next section for detail. 

The second set of sub critical values of g can be defined as the smallest value of g 
such that the function Y m+r \ vw acquires a double zero at u — 0. They are denoted as 
g r j™, and they are always greater than the corresponding critical values: g r j™ < g r j™. 

The locations of the critical values depend on the state under consideration, and 
can be determined approximately by using the asymptotic Y-functions discussed in 
the previous subsection. The values obtained this way are only approximate because 
for large enough values of g one should take into account the deviations of the Y- 
functions from their large J expressions. 

We will see in next sections that the critical values g r j™ play a crucial role in 
formulating excited states simplified TBA equations which take different form in 
each of the intervals g r j™ < g < g r j~ n ,m , r = 0, 1, . . . where g Q j™ = 0. The second set 
of g r j™ is not important for the simplified equations. The canonical TBA equations 
however require both sets because they take different form in each of the intervals 

r.m , , — r,m —r,m . . r+l.m n 1 

9 j,n <9< gin ; 9j,n <9<9j, n > r = 0, 1, . . .. 

Strictly speaking the integration contour in TBA equations also depends on g 
and the state under consideration. Nevertheless it appears that in simplified TBA 
equations the contour can be chosen to be the same for all values of g, and even for 
all two-particle states from the sl(2) sector if one allows its dynamical deformation. 
That means that with increasing g the contour should be deformed in such a way 
that it would not hit any singularity. This also shows that one should not expect 
any kind of non-analyticity in the energy of a state at a critical value of g. What 
may happen is that the critical values are the inflection points of the energy. 

Critical values of g for the Konishi state 

In this subsection we discuss in detail the critical values for Konishi state. To analyze 
the dependence of Y-functions on g one should first solve the BY equations with 
J = 2, n = 1, and then plug the u/s obtained into the expressions for Y-functions 
from appendix 

Solving the equations 

Y r+llvw {±~,g) = 0, r = l,2,..., (3.17) 

g 

we find that there are 7 critical values of g for g < 100 

g £f = {4.429, 11.512, 21.632, 34.857, 51.204, 70.680, 93.290} . (3.18) 
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Figure 3: On the left and right pictures Yi\ vw , Y 2 \ vw and Y s \ vw are plotted for the Konishi 



state at g$ ~ 4.5 and g c V ~ 11.5, respectively. Y2\ vw touches the n-axis at g = gcr , and 
has two real zeros for g^) < g < gcV ■ 



,(2) 



(1) 



Note that the distance between the critical values increases with g. The first critical 
value is distinguished because only Y 2 \ vw (±i/ g, g) vanishes there. For all the other 
critical values the function Y r _i\ vw (±i/ g, g) also is equal to zero 



Y , (+ l n {r ' lh 
9 







(r,lh 



Yr— i| wu ,(i , (72,1 
9 



, for r 



2,3, 



Then, solving the equations 

Y r+1 \ vw (0,g) = 0, r = 1,2,... , 
one finds the following 7 subcritical values of g for g < 100 

{4.495, 11.536, 21.644, 34.864, 51.209, 70.684, 93.292} . 



-(r,l) 
92,1 



(3.19) 



(3.20) 



(3.21) 



Note that the distance between a critical value and a corresponding subcritical one 
decreases with g. Again, at the first subcritical value only Y2\ vw (0,g) vanishes. For 
all the other subcritical values the function Y r _i\ vw (0, g) also acquires an extra double 
zero 



Y, 



r+l\vw 



(0, Ki 







Y, 



r—l\vw 



(r,lh 



, for r = 2, 3, 



(3.22) 



Once g crosses a subcritical value yjji the corresponding double zeros at u = split, 
and each of the functions Y r _i\ vw (0, g) and K+i|«w(0, g) acquires two symmetrically 
located zeros. As a result at infinite g all the YM|i>u>-functions have four real zeros. 
Moreover, one can also see that if g is between two subcritical values then all these 
zeros for all the functions are inside the interval [—2, 2] and approach ±2 as g — > oo. 
In Figures 3-5, we show several plots of YMi^-functions for the Konishi state. 

In what follows a Konishi-like state refers to any two-particle state for which only 
Yil^-function has two real zeros and all the other Y^^-functions have no zeros in 
the physical strip at small coupling. 
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Figure 5: On the left picture and right pictures Y 2 \ 



for the Konishi state at g, 

,(3) 



(3) 



21.6 and gif 



vwi ^3\vwj Y&\vw an d y^\vw are plotted 
34.9, respectively. Y 2 \ vw and Ym vw touch the 



u-axis at g = g cr • , and Y 3 \ vw and Y$\ vw touch it at g 



9 > 9cr ■ 



9cr ■ 



Y 2 \ vw has four real zeros for 



Critical values of g for some states 

Here we analyze the (/-dependence of Y-functions for several other states. 

We begin with the state with J = 5 and n — 1. This is the state with the lowest 
value of J such that both Yi\ vw and Y 2 \ vw have two zeros in the physical strip at small 
g. The critical and subcritical values are determined by the equations 

Y r+2lvw (±-,g)=0, Y r+2lvw (0,g) = 0, r = l,2,..., (3.23) 
9 

and we find the following values for g < 100 

g r 5 j = {6.707, 15.458, 27.233, 42.107, 60.101, 81.222} (3.24) 
gl* = {6.764,15.479,27.244,42.114,60.105,81.225}. 

Next, we consider the state with J = 8 and n = 1. This is the state with the 
lowest value of J such that Yu vw has four zeros, Y 2 \ vw has two real zeros and Y 3 \ vw 
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has two imaginary zeros in the physical strip at small g. Therefore, the critical and 
subcritical values are determined by the equations 

Y r+3lvw (± l -,g)=0, Y r+2lvw (0,g) = 0, r = l,2,.... (3.25) 

We find the following 6 critical and 7 subcritical values of g for g < 100 

g r 8 j = { -,9.157,19.561,32.985,49.505,69.143,91.909} (3.26) 
= {0.116,9.207,19.580,32.995,49.511,69.148,91.912}. 

1 3 

The reason why g 81 is so small is that the two imaginary roots of Y 3 \ vw that are in 
the physical strip at g — reach the real line very quickly. 

One might think that the first critical value increases with J. It is not so as one 
can see on the example of the state with J = 9 and n — 1. This is again the state 
such that Yi\ vw has four zeros, and Y 2 \ vw and Y 3 \ vw have two zeros at small g, and it 
has the following 6 critical values of g for g < 100 

g r g j = {6.970, 16.982, 29.935, 45.968, 65.114, 87.384} (3.27) 
g r 9 j = {7.052,17.006,29.947,45.976,65.119,87.388}. 

4. TBA equations for Konishi-like states 

As was discussed above to formulate excited state TBA equations one should choose 
an integration contour, take it back to the real line of the mirror plane, and then 
check that the resulting TBA equations are solved by the large L expressions for Y- 
functions. We begin our analysis with the simplest case of a Konishi-like state which 
appears however to be quite general and allows one to understand the structure of 
the TBA equations for any two-particle sl(2) sector state. To simplify the notations 
we denote the critical values of the state under consideration as giP. 

4.1 Excited states TBA equations: g < g^r 

Let us stress that the equations below are valid only for physical states satisfying the 
level-matching condition. Since some terms in the equations below have the same 
form for any N we keep an explicit dependence on N in some of the formulae. 

Integration contour 

The integration contour for all Y-functions but Y± is chosen in such a way that it 
lies a little bit above the interval Ke(z) G (— lm(z) = ^ in the middle of the 
mirror theory region, and penetrates to the string theory region in the small vicinity 
of z = u 2 /2 (the centre point of the mirror region). Then, it goes along the both 
sides of the vertical line to the string region real line, and encloses all the points 
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Figure 6: On the left picture the integration contour on the z-torus is shown. On the 
right picture a part of the integration contour on the mirror u-plane going from ±00 a little 
bit above the real line to the origin, and down to the string region is plotted. The green 
curves correspond to the horizontal lines on the z-torus just above the real line of the string 
region. The red semi-lines are the cuts of the mirror u-plane, and on the z-torus they are 
mapped to the part of the boundary of the mirror region that lies in the string region. 

z*k so that they lie between the mirror theory line and the integration contour, see 
Figure 6. In the mirror w-plane the contour lies above the real line, then it goes down 
at u — 0, reaches a minimum value and turns back to the real line. In the equations 
involving the functions Yq, Yq\ vw and Yq\ w it crosses the cuts of the mirror w-plane 
with Im(tt) = —gi and enters another sheet, see Figure 2. It is worth mentioning 
that the contour does not cross any additional cuts ^-functions have on the z-torus. 

Then, one uses the TBA equations for the ground state energy, and taking the 
integration contour back to the mirror region interval Im(z) = picks up N extra 
contributions of the form — log S(z*, z) from any term log(l + Yx)*K, where S(w, z) 
is the S-matrix corresponding to the kernel K: K(w,z) = t^I^ log S(w, z). In 
addition, one also gets contributions of the form — log S(w,z) from the imaginary 
zeros of 1 + Ym\ V w located below the real line of the mirror w-plane, see ( |3.16| ). 

Finally, the integration contour for Y±-functions should be deformed so that the 
points = Uk — - of the mirror u-plane lie between the interval [—2, 2] of the mirror 
theory line and the contour. Then, the terms of the form log(l — Y + )&K would 
produce extra contributions of the form + log S(u^, z) because Y + {u~^) = 00. In 
fact, this is important only if one uses the simplified TBA equations because in the 
canonical TBA equations Y±-functions appear only in the combination 1 — y-. Note 
also that Y±-functions analytically continued to the whole mirror w-plane have a cut 
(—00, —2] U [2, 00), and they should satisfy the following important equality which, 
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as was shown in [ 



"221 



], is necessary for the fulfillment of the Y-system 



Y + (u ± iO) = Y^(u T iO) for u G (—00, —2] U [2, 00) . (4.1) 

This equality shows that one can glue the two -u-planes along the cuts, and then 
l±-functions can be thought of as two branches of one analytic function defined on 
the resulting surface (with extra cuts in fact). We will see that the equality (|4.1| ) 
indeed follows from the TBA equations. 

Simplified TBA equations 

Using this procedure and the simplified TBA equations for the ground state derived 



in [22| |23|, one gets the following set of integral equations for Konishi-like states and 



g < g£> 

• M|w-strings: M > 1 , Y 0[w = 

\ogY M \ w = log(l + Y M -i\v,)0- + Y M+1 \ W ) -ks + 5 M i log _ j_ * s . (4.2) 

These equations coincide with the ground state ones. 

• iW>w-strings: M > 1 , Y \ vw = 

N 

\ogY M \ vw (v) = -SMiy^ogSju' - v) - log(l + Y M +i) *s (4.3) 

i=i 

1 - y_ 

+ log(l + F M -i|^)(l + Y M+1 \ VW ) -ks + 5 M i log — * s , 

1 — y + 

where uj = Uj — 1 , and the kernel s and the corresponding S-matrix S are defined 
in appendix |8.1j. For M = 1 the first term is due to our choice of the integration 



contour for l±-functions, and the pole of Y + at u = uj . 

It is worth pointing out that there is no extra contribution in ( |4.3D from any 
term of the form log(l + Yu vw ) * s. In general such a term leads to a contribution 
equal to — log S(ri — v)S(r 2 — v) where r\ and r 2 are the two zeros of 1 + Y k \ vw with 
negative imaginary parts. To explain this, we notice that if the zeros r^ k+1 ^ of Y k+ i\ vw 
lie outside the physical strip, then rj are r± = r[ k+1 ^ — - g and r 2 = r± + |, where 
Im(r[ fc+1 ^) < — K Since S(r — ^)S(r + |) = 1 the term log(l + Y k \ vw ) * s does not 
contribute, see Figure 7. 

On the other hand, if the zeros rj h+1 ^ of Y k+ i\ vw lie inside the physical strip then 
rj are related to rj fc as rj = rj k+1 ^ — ^, and the term log(l + Y k \ vw ) * s leads to the 
extra contribution equal to — £\ log S(rj k+1 ^ — v) where r^ 1 '' = rj k+1 ^ — |, see 
Figure 7. 
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Zeros of Y k+1 \ vw for g < g c 



Zeros of 1 + Y k u 



CM) 



log(l + Y k \ vw ) * s 
Zeros of Y k+1 \ vw for g cr < g < g cr 



iog(i + y fc | i;u; )*s 

Zeros of F fe+ i for <7 cr < 5 



CM) 



- r 
iX/\ 


r 




f 1 — ' 



log(l + lfc| w )*S 





cancel 



S(- r -±-u)S(r-±-u)^l 




do not 
cancel 



Figure 7: In the upper pictures the positions of zeros of Y k+ i\ vw and of 1 + Y k \ vw are 
shown for g < g cr . The contributions of log(l + Y k \ vw ) * s to the TBA equations cancel out 
for this case. The two pictures in the middle correspond to the situation when two zeros of 
Y k+ i\ vw enter the physical region |Imn| < | (depicted in yellow). Finally, the two pictures 
at the bottom are drawn for the case when two zeros of Y k+ ±\ vw are on the real line which 
corresponds to g cr < g cr < g. When g > g cr the contributions of log(l + Yu vw ) * s do not 
cancel anymore and lead to modification of the corresponding TBA equations. 
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We conclude therefore that at weak coupling and only for Konishi-like states no 
term log(l +Yk\ vw )*s gives an extra contribution to the TBA equations for t>w-strings. 

Let us also mention that the poles of S(uJ — v) cancel the zeros of Yi\ vw {y) at 
v — Uj, and ( |4.3|) is compatible with the reality condition for Y-functions. 



y-particles 4 

Y N 

log -±{v) = - log S Uy (u v v) + log(l + Y Q ) *K Qy , (4.4) 



3=1 



\ogY + Y4v) = -Y^log^-^*s( Uj ,v) (4.5) 



+ 2 lQ g 1 1 v™ * « ~ ^8 (1 + Y Q) *K q + 2 log(l + Y Q ) * K% * . 

1 + *l\w 



where we use the following notation 



* s(u v )= dt i g ^ v J ' I 



Then, S\ ty (uj,v) = Si y (z*j,v) is a shorthand notation for the S-matrix with the 
first and second arguments in the string and mirror regions, respectively. The same 
convention is used for other S-matrices. Both arguments of the kernels in these 
formulae are in the mirror region. 

Taking into account that under the analytic continuation through the cut \v \ > 2 
the S-matrix Si tV and the kernel Kq v transforms as Si tV — > l/Si ty and Kq y — > —Kq y , 
one gets that the functions Y± are indeed analytic continuations of each other and, 
therefore, the equality fl4.1|) does hold. 

It can be easily checked that the term on the first line in ( |4.5| ) is real, and 
this makes obvious that the equations for Y±-functions are also compatible with the 
reality of Y-functions. The origin of this term can be readily understood if one uses 
the following identity 



- ]og^p-*s( Uj ,v) = log S 1 {u j - v) 2 log * .s(„,, r) ( 1.6) 

<->2 



which holds up to a multiple of 2ni, and since S^i^u^) has a zero at u = v the 
integration contour in the second term on the r.h.s. of ( |4.6| ) runs above the real 
line. Then, the term log Si (uj — v) comes from the term — log (1 + Y\) * Ki, and the 
second term come from 2 log(l + Yi) * Kfgj * s. 



4 The equation for Y+Y_ follows from eq.(4.14) and (4.26) of p2) , and the identity 
Kq v * Ki — K®* — Kq-i . Another identity ★ s — \Kq y + \Kq — 5q\s — Kq^ * s , where 



s(u) — s(u — is useful in deriving Y-system equations for Y±. 
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Eq.(4.5) is very useful for checking the TBA equations in the large J limit where 
one gets 

N (^ U1 ) 2 i _i_ y,, 

logr + F_ = - V togi^* S (^,i;)+2togii^* a . 



, - l\w 

3=1 



Q-particles for Q > 3 



1 + 



Q = 2-particle 



N 11 / 1 \ 

log Y 2 = V log S( Uj -v)- log(l + — )(1 + — ) * 8 + 2 log 1 + — - * s (4.8) 



3 

In fact by using the p. v. prescription, one gets 

2 



(l + yf 

logy 2 = log ^j— r^jfr^ *p-»- s («) 



which makes obvious the reality of Y-functions. 
• Q = 1-particle 



~ / i X 

log Yt = ^ log & * v) - L £ * s + log ( 1 - — j Y 2 * s 
- 21og (1 - J (l - ] y 2 *#*s + logYi*#i*s 



- log (1 + F Q ) * (2KJ + K Q + K Q _ 2 ) * s - log(l + Y 2 ) * s . (4.10) 



All the kernels appearing here are defined in appendix |S.1| , and we also assume that 
Kq = and K-\ = 0. The reality of this equation follows from the reality of 

S ss (u — -, v) 2 x s (u — -) — x s {v ) x s (u + L ) — x s (v) i i 

h ( \ ~ ( i\ 1 / i i\ 1 — o ss \U ,V)O ss [U-\- ,V), 



which appears if one uses the representation ( |3.24| ) for S^. Note that S±(u, v) is 
defined through the kernel K\(u,v) by the integral 

S 1 (u,v) = ex P (2m / du'K x {u',v)) = -f (4.11) 

V J~oo J S SS {U+-,V) 
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and it differs from the naive formula 



i 

S ms (u - -,v)S ms (u + -, v) = v) x s {v) 2 , (4.12) 



which one could write by using the expression (|8.16| ) for the kernel K\. 

We see that the reality of Y-functions is a trivial consequence of these equations. 
Moreover, in the large L limit the simplified TBA equations do not involve infinite 
sums at all. As a result, they can be easily checked numerically with an arbitrary 
precision. We have found that for Konishi-like states the integral equations are solved 
at the large L limit by the asymptotic Y-functions given in terms of transfer matrices 
if the length parameter L is related to the charge J carried by a string state as 

L = J + 2. 

We expect that for all TV-particle states from the s((2) sector the relation between 
length and charge is universal and given by L = J + 2. 

There is another form of the TBA equations for Q-particles which is obtained 
by combining the simplified and canonical TBA equations. We refer to this form as 
the hybrid one, and the equations can be written as follows 

• Hybrid equations for Q-particles 

TV 

logY Q (v) = -J2 (lo g 5jg(ti if t;) -2togS*A*L(«7.*)) 

- L S Q + log (1 + Y QI ) * + 2 s* K^ X ' Q ) (4.13) 
+ 2 log (l + Y x \ vv} ) * s*K yQ + 2 log (l + Yq-iiw,) * s 

- 2]og±^-*s*K% x + log^-^^ Q + log (1 - i-) (1 - i_) *KyQ , 
where K^ x = 0, and Y \ vw = 0, and we use the notation 

/oo 
dt log S(uj -t-iO)* Kl2 x {t + iO, v) . (4.14) 
-oo 

The first term on the first line of ( |4.13| ) comes from log (1 + Yq>) ★ K%£, and the 
second one from —2 log is * K^ x . Eq. ( |4. 13|) is derived in appendix |8^4] . 

The energy of the multiparticle state is obtained in the same way by taking the 
integration contour back to the real mirror momentum line, and is given by 

N „ oo 7~Q 

E{n k] {L) = ~ / dZ E bg (1 + 

k=l J Q=l 

N 



E^- f du J:^ d -L^(l + Y Q ), (4.15) 



2n du 

k=l •' Q=l 
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where 



8 k = igx (z* k ) - igx + (z* k ) - 1 = igx s (u k ) - igxj(u k ) - 1 , (4.16) 



is the energy of a fundamental particle in the string theory, see appendix |0] for 
definitions and conventions. 

For practical computations the analytic continuation from the mirror region to 
the string one reduces to the substitution x^iu) — » xf^^u) = x s {u ± -Q) in all the 
kernels and S-matrices. Then, as was discussed above the string theory spectrum 
is characterized by a set of N real numbers u k (or z* k ) satisfying the exact Bethe 
equations (|3.7|). We assume for definiteness that u k are ordered as U\ < ■ ■ ■ < u N . 

Finally, to derive exact Bethe equations one should analytically continue Y\ given 
by either eq. ( |4.1(J| ) or ( |4. . We find that it is simpler and easier to handle the exact 



Bethe equations derived from the hybrid equation Q4.13 ) for Y\. In the appendix [|]6] 



we also derive exact Bethe equations from the canonical equation for Y\. 
Exact Bethe equations 

Now we need to derive the integral form of the exact Bethe equations (|3.7|) . Let us 
note first of all that at large L eq. ([3.71) reduces to the BY equations for the sl(2)- 
sector by construction, and the integral form of (|3.7|) should be compatible with this 
requirement. 

To derive the exact Bethe equations, we take the logarithm of eq.(|3.7|), and 



analytically continue the variable z of Y\{z) in eq. (|4. 1 3|) to the point z* k . On the 



mirror w-plane it means that we go from the real w-line down below the line with 
Im('u) = — - without crossing any cut, then turn back, cross the cut with Im(-u) = — 
and |Re(ii)| > 2, and go back to the real w-line, see Figure 6. As a result, we should 
make the following replacements x(u — -) — > x s (u — -) = x(u — -), x(u + -) — > 
x s (u + 1 ) = \jxiu + -~) in the kernels appearing in ( [4. 131 ). 

The analytic continuation depends on the analytic properties of the kernels and 



Y-functions, and its detailed consideration can be found in appendix p. 51 As shown 
there, the resulting exact Bethe equations for a string theory state from the sl(2) 
sector can be cast into the following integral form 



N 



ni(2n k + 1) = log>i,(ufc) = iLp k - ^ lo S ^(2) ( u j > u k) ( 4 -17) 



2k x l _ 



N N 2 ' X J ~ ' 

+ 2 log Res(S) * Kltiuj, u k )-2Y^ log {uj -u k --) 

j=l j= i 9 •'•/,• 

+ log (1 + Y Q ) * (tfji. + 2 s * K^ 1 *) + 2 log (1 + Y Avw ) ★ (a* K yU + s) 

~ 21 °s j^^* s * Kl vt + log Y_7 "X*^ 1 + lo s ( x - Y^i 1 - y^)* k vu , 
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where we use the notations 



/+oo 

dtlog S(u--t)(t-u) Kl l : x {t,v), (4.18) 
-oo 

s{u)=s{vT). (4.19) 

The integration contours in the formulae above run a little bit above the Bethe roots 
u jj Pk = i£-Q(z*k) = —i log X . S | M ^|| is the momentum of the k-th particle, and the 



second argument in all the kernels in ( 4.17|) is equal to The first argument we 



integrate with respect to is the original one in the mirror region. 

Taking into account that the BY equations for the st(2) sector have the form 



N 

■Ki(2n k + 1) = i Jp k - ^ lo S S li\l) ( u j> u k) , 

j=i 



and that Yq is exponentially small at large J, we conclude that if the analytic con- 
tinuation has been done correctly then up to an integer multiple of 2ni the following 
identities between the asymptotic Y-functions should hold 



N _ N 2i X J- A 

TZ k = 2ip k + 2^2 logRes(S') * K^ x (Uj , u k ) - 2 ^ log (v,j - u k - — ) 

3=1 3=1 

I — y 

+2 log (1 + Y Mvw ) *(s* K yU + S) - 2 log — — * s * Kill 



k 



g Xj 



+ log l —^ r ±K 1 + bg (1 - -L) (1 - -L) *K vU = . (4.20) 

For N = 2 and u\ = —112 one gets one equation, and by using the expressions for the 
Y-functions from appendix |S.3| one can check numerically 5 that it does hold for any 
real value of u\ such that only Yi\ vw has two zeros inside the strip |Imu| < 1/g. 

4.2 Excited states TBA equations: g^r < g < gcr 

In this subsection we consider the TBA equations for values of g in the first critical 
region g^ < g < gcr ■ In this region in addition to the two real zeros of Y\\ vw at u 
two zeros of Y 2 \ vw enter the physical strip. We denote these zeros r,-. 



3 ' 



5 Since K^ x (u,v) has a pole at u = v with the residue equal to —5^7 the terms of the form 
2f * K^ x can be represented as 2/ * K^ x = 2/ * p ,„. R^wx + f( u k) which is useful for numerics. 

6 Canonical TBA equations take different forms depending on the reality of the zeros, and, 
therefore, one has to divide the region into two subregions: < g < gil^ and g^ < g < gc?, see 



the next subsection and appendix B.6 for detail 
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Simplified TBA equations 

We first notice that for g$ < g < gf) the function Y 2 \ vw has two zeros in the 
physical strip. Therefore, as was discussed above, the contribution to the simplified 
TBA equations coming from the zeros of Yi\ vvi and 1 + Y\\ vw does not vanish. If 
g < g^cr no deformation of the integration contour is needed because all these zeros 
are on the imaginary line of the mirror region. As g approaches g^/ the two zeros of 
Y-i\vw approach u = 0, and at g = g^ they both are at the origin. As g > gil^ the 
zeros become real, located symmetrically, and they push the integration contour to 
be a little bit below them. In addition to this at g = g$ the two zeros of 1 + Yi\ vw 
with the negative imaginary part reach the point u = —i/g, and as g > g$ they 
begin to move along the line Im(-u) = — 1/g in opposite directions. As a result, the 
integration contour should be deformed in such a way that the two zeros of 1 + Yi\ vw 
would not cross it. Thus, the zeros of 1+Yi\ vw always lie between the real line and the 
integration contour, and the terms of the form log(l + Yi\ vw ) * K produce the usual 
contribution once one takes the contour back to the real line. Let us also mention 
that the points rj of the mirror u-plane are mapped to the upper boundary of the 
string region on the z-torus. 

Using this integration contour, one gets the following set of simplified TBA 
equations for Konishi-like states and g^ < g < gc 



• M|ty-strings: their equations coincide with the ground state ones (4.2). 

• M|t>ty-strings: M > 1 , Y \ vw = 

N 2 

log Y M \ vw {v) = -5 mi ^ ^gS(uj -v)- 5 M2 ^2 1°S 5 ( r 7 - v) (4.21) 

i=i i=i 

1 - Y- 

+ log(l + Y M -x\ vw ){l + Y M +i\vw) *s + 6 M1 log — *s - log(l + Y M+1 ) * s . 

1 ~~ Y + 

The first term is due to the pole of Y + at u = uj , and the second term is due to the 
zeros of 1 + Yi\ vw at u = rj . 

• ^/-particles 

Y - 

\og^(v) = -J2^gS Uy ( Uj ,v) + \og(l + Y Q )^K Qy , (4.22) 
i=i 

\ogY + Y4v) = - log * s( Uj , v)-2J2 log S(rj - v) (4.23) 

3=1 2 3=1 

+ 2 l°g 1 1 t Y ' lVW *s-\og(l + Y Q )*K Q + 2\og(l + Y Q )*Kg l *s, 

1 + *l\w 
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where the second term on the second line is due to the zeros of 1 + Yi\ vw at u — r- 
• Q-particles for Q > 3 



log Yq = log - , v j- y «- 1| 7' j- - s . (4.24) 



In fact the p.v. prescription, see appendix |8.7| , is not really needed here because for 
Q = 3 the double zero of Y 2 cancels the zeros of Y 2 \ vw , and for Q > 4 everything is 
regular. Thus, the formula works no matter if the roots rj are real or imaginary. 

• Q = 2-particle 

'l + 



log Y 2 = -2 J2 S(rj -v) + log V ^ 5 (4.25) 

This makes obvious the reality of Y-functions because the double zero of Y 2 at v — rj 
is cancelled by the pole of S(rJ — v). 

Hybrid equations 

One can easily see that the simplified equation for Q = 1-particles is the same as 
eq.( |4.10|) in the weak coupling region g < . Thus, we will only discuss the hybrid 



equations for Q-particles. Strictly speaking, their form is sensitive to whether the 
zeros of Y 2 \ vw are complex or real, and the first critical region is divided into two 
subregions: g^ < g < g^ and g^ < g < gir . On the other hand, to derive the 
exact Bethe equations we only need the hybrid Q = 1 equation which, as we will see, 
takes the same form in both subregions, and, moreover, for any g > g^r . 

For g^ < g < g$ the function Y 2 \ vw has two complex conjugate zeros in the 
physical strip which, therefore, lie on the opposite sides of the real axis of the mirror u- 
plane. Thus, the zero T\ with the negative imaginary part lies between the integration 
contour and the real line of the mirror region. Taking the integration contour back to 
the real line produces an extra contribution from this zero. Since 1 + Yi\ vw (rj) = 0, 
and Kj_ = 0, the terms 2 log (l + Y^ vw ) * s* K yl and log (l — ^-) (l — p-) *K y i 
in the hybrid equation ( |4.13| ) lead to the appearance of — 2 log S * K yl (r J ,v) and 
2logS y Q(ri,v), 7 respectively. The integration contour for the * -convolution in the 
term — 2 log S*K y x(rJ, v) is the one for F±-functions and it should be taken to the 
mirror region too. Since the S-matrix S(r^ — t) has a pole at t = r 1; this produces 
the extra term — 2logS y Q(ri,v) that exactly cancels the previous contribution from 
log (l — p-) (l — p-) * Kyi. Thus, the only additional term that appears in the hybrid 

TBA equation (p^) for g£) < g < gS) is -2J2^ogSiK vl (rJ,v). 



''Note that the S-matrix S v q is normalized as S y Q(±2,v) = 1 
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As g > gir the two zeros of Y± (and Y 2 \ vw ) become real and located a lit- 
tle bit above the integration contour. Therefore, the only extra contribution to 
the TBA equation comes from the term 2 log (l + Yi\ vw ) * sicK y i and it is again is 
-2^ogS*K yl (r~,v). 

Thus, we see that in the first critical region gil < g < gc? the hybrid Q = 1 
equation takes the following form independent of whether g < gt/ or g > g$ 
• Hybrid Q = 1 equation for g$ < g < g$ 



\ogY 1 {v) = -J2 (kg -2logS*K£, x (uJ,vj) -2j2^gS*K yl (rJ 

3=1 3=1 

- L E x + log (1 + Y Q ) * (K% 2) +2s* tf*" 1 ' 1 ) + 2 log (1 + Y 1{vw ) *s*K yl 



-2tog[-^*«*i^ + bgi-^*Jfi + bg(l-^)(l-^)*lf I , 1 . 

(4.26) 

Note also that for g > g^ one can also use the p. v. prescription in the terms 
log S * K y i(rJ , v) and log (l — p-J (l — y-j *K y i because the extra terms ± log S y i 
cancel each other. 

Exact Bethe equations: g^ < g < gcr 

The exact Bethe equations are obtained by analytically continuing logYl in ( f4.2G| ) 
following the same route as for the small g case, and they take the following form 

N 

ni(2n k + 1) = log Y u (u k ) =iLp k -^ lo § S li(l) K'> u k) ( 4 -27) 

3=1 

- _ - 2% x 3~^ 
+ 2 1o § Res(5) -k K l ^ x ( Uj , u k ) - 2 Y log (uj - u k ) _ _ \ 

3=1 3=1 9 X i Xk 

2 

- 2 2j fog S*K yU (rJ,u k ) - log S(rj - u k )) 
3=1 



+ log (1 + Y Q ) * [K^ + 2s* S- 1 ' 1 - ) + 2 log (1 + Y llvw ) *{s* K yU + s) 
- 2 log * 8 * iC; + log * x + log (1 - i-) (1 - -L) * 



We recall that in the formulae above the integration contours run a little bit above 
the Bethe roots Uj, and below the dynamical roots rj. The exact Bethe equation 
also have the same form no matter whether g < g$ or g > g^ . 

We conclude again that the consistency with the BY equations requires fulfill- 
ment of identities similar to ( |4.20 ) that we have checked numerically. 
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Integral equations for the roots rj 

In the first critical region the exact Bethe equations should be supplemented by 
integral equations which determine the exact location of the roots r 3 -. Let us recall 

i 



that at rj = r j ± - the function Y\\ vw satisfies the relations 



Y 1]vw (rf) = -1. (4.28) 

These equations can be used to find the roots rj. They can be brought to an integral 
form by analytically continuing the simplified TBA equation ( |4.21| ) for Yi\ vw to the 
points r^. The analytic continuation is straightforward, and one gets for roots with 
non-positive imaginary parts 

- ' i _ y_ 



IlS(, 

k=l 



u k - Tj) exp log(l + Y 2 \ vw ) -k s + log — *s- log(l + Y 2 ) ★ s) = -1 ,(4.29) 



where we used the exponential form of the TBA equation, and took into account that 
Yz\vw(rj) = Y 2 (rj) = Y±(rj) = 0. One can easily check that if rj is imaginary, that 
is g < gir\ then the product of S-matrices in (|4.29f) is negative, while the exponent 
is real. As a result, taking the logarithm of ( L29|) does not produce any new mode 



number. If rj is real, that is g > §cr, then both factors in ( 4.29Q are on the unit 
circle. Their phases however are small (as they should be due to the continuity of 
the roots as functions of g), and one does not get any mode number either. 

Let us also mention that eq. ( |4.29| ) can be also used to find the critical values of 
g. If one sets rj = —i/g, one gets an integral equation for gil\ and the case of rj = 



gives the one for gil\ The equation ( 4.29|) is better for numerical computation than 



( |4.28 ), because all the Y-functions are evaluated on the mirror real axis. 
4.3 Excited states TBA equations: g^ < g < ^c™ 



lev 



In this subsection we propose the TBA equations for Konishi-like states for the gen- 
eral case where g^ < g < gi™^ and m > 2. Then, the functions Y x \ vw , . . . , Y m _i\ vw 
have four zeros, and Y m \ vw and Y m+ i\ vw have two zeros in the physical strip. The only 
two roots that can be imaginary if g < g^f* are A m . As was discussed in section 



3, these zeros can be written in the form (P -14 ) 



{Uj , rf} , {rf , rf } , . . . , {rf " 1} , rf , {,>>} , {rf , (4.30) 

where we indicate only those zeros which are in the physical strip and important for 
formulating the TBA equations. We also recall that the functions below have zeros 

(k) 

at locations related to 

Y ± (rf)=0, l + Y k \ vw {rf +1) ± l -)=Q 1 Y k+l (rf = , k = 1 , . . . , m . (4.31) 
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Simplified TBA equations: g cr ' < g < g, 



cr 



The necessary modification of the integration contour is obvious, and simplified equa- 
tions take the following form 

• M|w;-strings: their equations coincide with the ground state ones (4~2|). 

• M|fw-strings: M > 1 , Y \ vw = 

2 

\ogY M]vw (v) = ~J2 [log^r^ - v ) + log^(rj M+2) '- - v)] (4.32) 



j=i 



1 - Yl 



+ log(l + Y M -i\vw)(l + Y M +i\vw) *s + 5 M i log — * s - log(l + Y M +i) ★ s 

i — y _i_ 



where we identify = Wj, and assume that log S(rj k ^' — v) is absent in the sum 
on the first line if k > m + 2. For M = 1 the term log S^Wj — v ) is due to the pole of 

— (3) — 

Y + at m = it • , and the second term is due to the zero of 1 + Y 2 \ vw at u = r y - '' . For 
M > 2 both terms on the first line are due to the zeros of 1 + Y k \ vw at u = r)- k+1 '~ 



for k = 1, . . . , to. 



y-particles: their equations coincide with eqs. ( [4.22] ) and ( [4.23| ) for the first critical 



region gil> < g < gi? . 



Q-particles for Q > 3 

2 ( I + 



log Y Q = -2J2 log S(rf - «) + log a . (4.33) 

j=i { + yg_Jv A + Y Q+1 ) 



The p.v. prescription is again not really needed because the four zeros 

of Yq_i\ vw are cancelled by the double zeros of and Vq+i, and the formula works 

no matter if the roots A m are real or imaginary. 



The equations below coincide with the corresponding eqs. ( |4. 25| ) , ( [4 . 1 0| ) , ( |4. 26| ) and 



( |4.27| ) for the first critical region g^) < g < gi) with the replacement Tj i— . 

• Equations for Q = 2 and Q = 1-particles. 

• Hybrid Q = 1 equation. 

• Exact Bethe equations. 

According to these results, only some of the simplified equations for Q-particles 
and t>u>-strings change their form when crossing the m-th critical point (to > 1). 
However, since all Y-functions are coupled with each other, the existence of higher 
critical points still affects the equations above but in a less direct way. 
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Figure 8: In the left and right pictures the profiles of the asymptotic Y\ and Y 2 functions 
for the Konishi state at g = gil^ = 4.429 are depicted, respectively. 

Integral equations for the roots r^: < g < g^ +V> 

In this region the exact Bethe equations should be supplemented by integral equations 
which determine the exact location of the roots rj , k — 2, . . . , m + 1. The integral 
equations are obtained by analytically continuing the simplified TBA equations ( |4.32| ) 
for Yk\vw to the points rj fc+1 ^ ± , and using the conditions 

Y klvw (r? +1)± ) = -1 ■ (4.34) 

The analytic continuation is again straightforward, and one gets for roots with non- 
positive imaginary parts the following equations 



rf : J] S{u n - rf )S(r® - r f) x (4.35) 

n=l 



xoxp( log(l + Y 2 \ vw ) ★g + log ^y *s- log(l + Y 2 ) * s ) = -I. 



2 



n ■ - : 11 ^ - r? +1) )S(r£ + V - r f x (4.36) 



.(*+!) . JJ (<? ^(ft) _ r (k+l)- ) q (r (k+2) _ r (fc+l)^ 

n=l 

x exp (log(l + + ife+i|im;) * 5 - log(l + Y k+1 ) *s) = -1, 

The continuity condition of the roots as functions of g again guarantees that one 
does not get any new mode numbers. 

4.4 Further evidence for critical values of g 

As was mentioned in the previous section the exact location of the critical values can 
be found only by solving the TBA and exact Bethe equations. One may wonder if the 
critical values can be absent at all due to the contribution of Yq functions. Indeed, 
lQ-functions at g ~ g cr are not small as one can see from Figure 8. However, their 
contribution to the quantities of interest is still small. To illustrate this, in Figure 9 
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Figure 9: In the left picture the profiles of the asymptotic Y^ vw (blue) and Yi\ vw after 
the first iteration (purple) are depicted. In the middle and right pictures similar profiles 
for Y + and Y_ are presented. 

we present the plots of asymptotic Y±, Yi\ vw and Y±, Yi\ vw (only for u > because 
all Y-functions are even) obtained after the first iteration of the simplified TBA 
equations by taking into account the contribution of the first eight Y^-functions for 
g = 4.38671, A = 759.691. The blue curves are plots of the asymptotic Y-functions. 
As we discussed in this section, at the first critical value Y^ vw (0) = — 1, and at the 
first subcritical value Y±(0) = 0. One can see from the plots that the influence of 
YQ-functions at u ~ is extremely small 8 . This is of course expected because the 
finite contribution of YQ-functions cannot balance the infinities originating from the 
zeroes of Y-functions. We also see from the plot of Y\\ vw that there is a tendency for 
the actual critical value to be higher than the asymptotic one. 

To show that the zeroes of Y-functions unavoidably enter the physical strip when 
the coupling increases, we also provide a plot of asymptotic Yg-functions for large 
value of A = 10 8 , see Figure 10. One can see that the functions are very small 
for almost all values of u except those close to ±2. Thus, at strong coupling all 
Y-functions are well-approximated around u = by their asymptotic value. On the 
other hand, if the zeroes are outside the physical strip then they are purely imaginary, 
and their positions can be determined by analytically continuing the corresponding 



TBA equations, see e.g. eq. (|4.29|) . Assuming the roots rj are imaginary, at large g 



the main contribution of the TBA kernel s appearing in ( f4.29|) originates from the 



region around u = 0. In this region, as was explained above, it is legitimate to use 
the asymptotic solution which implies that the roots rj are real and close to ±2. 
This leads to a contradiction with the assumption that the roots rj are imaginary at 
sufficiently large g. Thus, one concludes that for large g the roots are real and close 
to ±2. 

5. TBA equations for arbitrary two-particle 5 ((2) states 

Our consideration of the TBA equations for Konishi-like states can be easily used to 



8 Note, that for g far enough from the critical value, e.g. g = 2, the contribution of YQ-functions 
around u = results in a more visible change. 
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Figure 10: Profiles of the asymptotic Yq, Q = 1, 2 are depicted for A = 10 8 . 



formulate TBA equations for an arbitrary two-particle state from the sl(2) sector. 
One starts again with analyzing the structure of zeros of yM|i™-functions at weak 
coupling. Then, as was discussed in section 3, for every state (J, n) there is a number 
m such that the first m — 1 Y^^-functions have four zeros, both Y m \ vw and Y m+ i\ vw 
have two zeros, and all Y^\ vw -i unctions with k > m + 2 have no zeros in the physical 
strip. This is exactly the structure of zeros of Y^w-functions for a Konishi-like state 
with the coupling constant being in the m-th critical region: < g < g^ 1 . 

Thus, at weak coupling the simplified TBA equations for this (J,n) state have 
exactly the same form as the TBA equations for a Konishi-like state in the m-th 
critical region. As to the canonical TBA and exact Bethe equations, one should take 
into account that real zeros r!j of Yfci^-functions are located outside the interval 



[—2, 2] for small g for most of the s[(2) states. In particular, if the zeros of Y" 2 



(2) — 

are outside the interval then the corresponding zeros rj of 1 + Yi\ vw 
the cuts of Sl® x (r^ ,v), and since the integration contour runs below these zeros, 
one should add +i0 to them in all the expressions. In addition, the integration 
contour for Kj_-functions runs over the interval [—2, 2] (and a little bit above it), and 



j -» 2\vw 

are located on 



(2) 

if \rj | > 2 then there is no singularity of the integrand and the p. v. prescription 
is not necessary. Strictly speaking, this means that in addition to the critical and 
subcritical values of g we discussed above, one should also consider such values of g 
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that the zeros r - of Y 2 \ vw are equal to ±2, and distinguish the cases with the zeros 
being inside and outside the interval [—2,2]. It is straightforward to do it, and we 
refrain from presenting explicit formulae for these cases here. 

Increasing g, one reaches the first critical value of the (J,n) state, and 

the TBA equations would have to be modified. It is in fact clear that in the region 
^r,m+i < g < g r J { n ,m+ the TBA equations for the (J,n) state coincide with the ones 
for a Konishi-like state in the m + r-th critical region. 



6. Remarks on the Y-system 

For reader's convenience, in this section we summarize the most essential properties 
of the Y-system implied by the TBA equations under study. 

The kernel s has the uniquely defined inverse s _1 which acts as the following 
operator 

(/ * s- 1 )^) = lim [f(u + --ie) + f(u-- + ie)] . 

e->0+ Q 9 

Applying s _1 to the set of the simplified TBA equations, one can easily derive the 
Y-system equations [pT] for — 2 < u < 2. Since the Y-functions appear to be non- 
analytic, for \u\ > 2 (here and in what follows we mean u G (— oo, —2) U (2, oo) by 
\u\ > 2) one gets instead new equations which, as we will explain below, encode 
the jump discontinuities of the Y-functions across the cuts p2| |. To exemplify this 
statement, it is enough to consider Y\\ w and Yi-functions. 

We start with eq. (|4.2j) for Yi\ w . Applying to this equation, one finds 



1 i 

*}$ {«+*-- ( u - \ + <0 ) = i 1 + r^ (u) ' H < 2 (6 - 1} 

9 9 \ / i (Q) 

Y$ (u+- g - <0)yW (u - I + zO) = 1 + , \u\ > 2 . (6.2) 

We stress that these equations are unambiguously derived from the TBA equation 
and the e-prescription on the left hand side is fixed by that of We recall that in 
the TBA equations the functions Y± have their support on [—2,2]. 

The TBA equation ( [4. 2D for Y\\ w shows that this function has branch points 
located at u — ±2 ± K Since we are dealing with the mirror w-plane, it is natural 
to choose the cuts to run from ±oo to ±2 ± - parallel to the real axis. In the same 
way the TBA equations show that various Y-functions have branch points located 
at ±2 ± 1 Q, Q = 0, 1, 2, . . . , oo, and, therefore, all the cuts of all the Y-functions can 
be chosen to be outside the strip |Rew| < 2, and running parallel to the real axis. 
Then, Y-functions are analytic in the strip |Rew| < 2, in eq.( pM| ) the e-prescription 
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can be dropped, and u can be considered as a complex variable taking values in the 
strip 

1 i 

(« + (« - ~) = (l + *$) («) > | Re m | < 2 . (6.3) 



Thus, we see that with this (and only this) choice of the cuts the Y-system takes its 
standard form. 

In fact, analytically continuing Y-functions outside the strip |Reu| < 2, one 
concludes that eq.( |6.3| ) is valid for all complex values of u but those which belong to 
the cuts. Approaching then the real axis for \u\ > 2, say, from above, one arrives at 
the following prescription [2{| 9 



1 - 



i 



Y$ + iO)Y$ (u- l -+iO) = (l + Y^)-^{u + ,0) , \u\ > 2 . (6.4) 

"+ 

Such an analytically continued equation can be used to determine the jump discon- 
tinuities of the Y-functions across the cut. From the compatibility of the equation 
(gg) with eq.flOD, we find that 



*£ } («+S + <°) l ~y 



(a) 



and analogously 

In V o / 



(u + iO), |w|>2, (6.5) 



— (u-iO), |u|>2. (6.6) 



riot) 



n 

We also find that these discontinuity equations are consistent with the relation (|4.1| ) 
which follows from the TBA equations for Y±. 

Thus, the TBA equations tell us that, because of the absence of analyticity, the 
Y-system equations must be supplemented by proper jump discontinuity conditions. 
This was one of the important observations made in 

Finally, we recall that for Yi one finds 

1 



yd) J V yV 



Ki(u + --z0)ri(u--+i0) = - — — , |u|<2 

9 9 l+yr ( 6 .7) 



yi(u + --zO)Ki(u--+iO) = r , \u\ > 2 



9 9 ! + i 



9 Naively one could try to define a new "inverse" operator s^ 1 , 

(/*s (3 1 )(m) = lim e ^ + + ^ + ze) + i + ie)] > and claim that applying it to the simplified 
TBA equation, one gets eq.( |6.4[ ) for all values of u. The problem with this is that the operator s^ 1 



is not inverse to s. In fact it annihilates s: s-k Sq 1 = 0. 
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where the explicit form of the quantity A is given in [22l |2^| for the ground state 



case, and in the excited state case it can be extracted from eq.( |4.10|) or obtained from 



the ground state one by using the contour deformation trick. Once again, the second 
equation here determines the jump discontinuity of Y"i across the cut. However, a 
new feature here is that the jump discontinuity, that is A, does depend on a state 
under consideration! 

To summarize, the Y-system exhibits the following properties dictated by the 
underlying TBA equations 

• The Y-system is not analytic on the w-plane and for this reason it must be 
supplemented by jump discontinuity conditions; 

• In general, jump discontinuities depend on a state of interest and they can be 
only fixed by the TBA equations; 

• Different Y-functions have different cut structure. In total there are infinitely 
many cuts on the mirror theory w-plane with the branch points located at 
±2 ± 1 Q, Q = 0,1,2,..., oo. As a result, the Y-system lives on a Riemann 
surface of infinite genus. 



Needless to say, these intricate analyticity properties discovered in |22|, £9| render 
the AdS/CFT Y-system rather different from its known relativistic cousins and, for 
this reason, make it much harder to solve. 



7. Conclusions 

In this paper we have analyzed the TBA equations for the AdSs x S 5 mirror model. 
We provided an evidence that for any excited string state and, therefore, for any 
M = 4 SYM operator there could be infinitely many critical values of 't Hooft's 
coupling constant at which the TBA equations have to be modified. It is a demanding 
but also a rather challenging problem to locate the exact position of the critical points 
and it is similar in spirit to determination of the exact positions of Bethe roots. At 
the same time, we also want to give a word of caution. Our approach is based on the 
optimistic assumption that the analytic structure of the exact TBA solution emulates 
the one of the large L asymptotic solution. The possibility that the exact solution 
might develop new singularities in comparison to the asymptotic one would lead to 
even more complicated scenario than we described here. 

One could also speculate about the physical origin of critical points. It is known 
that the sl(2) sector is closed to any order of perturbation theory. One possibility 
would be that the first critical value of g of an 5 [(2) state is the one where this state 
begins to mix with states from other sectors of the theory. It would be interesting 
to understand if this is indeed the case. 
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The TBA equations and the contour deformation trick we have formulated allow 
one to discuss many interesting problems, and we list some of them below. 

Concerning the issue of critical points, it would be very interesting to analyze 
the TBA equations, and compute numerically the scaling dimension of the Konishi 
operator in the vicinity of and beyond a critical point. The simplified TBA equations 
seem to be much better suited for such an analysis than the canonical ones. 

Then, one should solve analytically the TBA equations for two-particle states 
at large g. The large g expansion should contain no logg which follow from the BY 
equations [^, It should also fix the coefficient of the subleading 1/g term in the 
energy expansion. There are currently two different predictions for the coefficient 
[j32| |3~I| obtained by using some string theory methods and relaying on certain as- 
sumptions. Thus, it is important to perform a rigorous string theory computation of 
this coefficient. 

Recently, the TBA approach has been applied to obtain the 5-loop anomalous 
dimension of the Konishi operator by a combination of analytical and numerical 
means [7D| and the corresponding result was found to be in a perfect agreement with 



the one based on the generalization of Liischer's formulae |T5[ . This constitutes an 
important test of the TBA equations we propose in this paper (the hybrid equations). 
It would be nice to support this numerical agreement by an analytic proof. 



The TBA equations we have proposed are not valid for the two-particle (J, ^jp) 
state. In the semi-classical string limit (7—7-00 and J/g fixed it should correspond to 
the folded string rotating in S 2 [53~|. It would be interesting to analyze these states 
along the lines of our paper, write TBA equations, and solve them. 

We have discussed only two-particle states. It is certainly important to generalize 
our analysis to arbitrary iV-particle 5 1(2) states. For J = 2 the lowest energy N- 
particle state is dual to the twist-two operator that plays an important role in field 



theory, see e.g. [64 



It would be also of interest to consider the one-particle case at large g and finite 
J/g. Its energy should match the string theory result |J5] . 

Let us finally mention that one should also consider other sectors and exhibit 
the psu(2,2|4) invariance of the spectrum. 
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8. Appendices 

8.1 Kinematical variables, kernels and S-matrices 

All kernels and S-matrices we are using are expressed in terms of the function x(u) 

1 



x(u) = -(u — i\/4 - u 2 ), Imi(n) <0, (8.1) 

which maps the u-plane with the cuts [—00, —2] U [2, 00] onto the physical region of 
the mirror theory, and the function x s (u) 



Xs (u) = l{l + \l-\), \x a (u)\>l, (8.2) 



2 V" V " v? 

which maps the w-plane with the cut [—2, 2] onto the physical region of the string 
theory. 

The momentum and the energy Sq of a mirror Q-particle are expressed in 
terms of x(u) as follows 

_ i i ~ x ( u ~ ~Q) 

Pq = gx(u- -Q) - gx(u + -Q) + iQ , £ Q = log . f— > . (8.3) 
9 9 x\u ~r ) 

The kernels act from the right, and the three types of star operations used in 
this paper are defined as follows 

/oo p2 
duf(u)K(u,v) , f*K(v)= duf(u)K(u,v), 
00 J -2 

f*K(v)=y + jf°) duf(u)K(u,v). (8.4) 

The TBA equations discussed in this paper involve convolutions with a number of 
kernels which we specify below, see also |22 for more details. First, the following 
universal kernels appear in the TBA equations 



1 d 1 gQ u ~ l ~ 

^lTu logSQ{u) = -,W^' ^ = ^' 



9 

1 d M ~ l 
K M n(u) = — — log S M n{u) = K M+N (u) + K N _ M (u) + 2 ^ k n~m+2j{u) , 

711 u j=i 

M-l 

Smn(u) = Sm+n{u)Sn-m(u) Y\ S N - M +2j(u) 2 = S NM (u) . (8.5) 
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Then, the kernels K± y are related to the scattering matrices S± y of Q- and y±- 
particles in the usual way 



K Qy (u,v) 
K® y (u,v) 



1 d 

2ni du 



1 d 

2ni du 



logS® y (u,v), S® y (u,v) 



log S?(u,v), Si y (u,v) 



x(u — i®) — x(v) 



x(u + ij) 



x(u + i&) — x(v) \ x(u — i^) ' 

x(u — l — ) 

v 9 > x[v) 



X[U 



x(u + i—) 



These kernels can be expressed in terms of the kernel Kq, and the kernel 



1 d 



1 y/A^ 



K(u,v) = — log S(u,v) = 

v 1 2m du v ' 2ni ^A-u 2 u-v 



, S(u,v) = 



x{u) — x(v) 
x(u) — l/x(v) 



,(8.6) 



as follows 



= l(K Q (u-v)±K Qy (u,v)y 



(8.7) 



where Kn y is given by 



K Qy (u,v) 
SQy(u,v) 



log S Qy (u, v) = K(u - -Q, v) - K(u + -Q, v) , (8.8) 

S Qy (u, v) _ x(u - |g) - xjv) xju + |g) - ^ _ Sju - |g, tQ 
S?(«, «) x(« - jQ) - ifo *(« + JQ) - x(v) S(u + jQ, W ) • 



The S-matrices 5^ and Sq^ (and kernels and Kq v ) can be easily continued to 
the string region by using the substitution x(u±i®) — > x s (u±i®), and the resulting 
S-matrices are denoted as S±* y and Sq, v . Notice that Sq^ v (— oo, v ) = 1. One can 
also replace x(v) by x s (v) in the formulae above. Then, one gets the S-matrices and 
kernels which are denoted as Kq^ (ms for mirror-string) and so on. 

The following kernels and S-matrices are similar to K± y 



K yQ (u,v) = ^--^-logS yQ {u,v), S yQ {u,v) 



K yQ (u,v) 



1 d 

2ni du 



log Sf(u,v), Sf( 



u,v) 



x(u) — x{v + i—) 



x(u) — x{v — ij) \ x(v + ij) ' 



x(v-i%) 



i 

x(u) 



x(v — i—) 

v 9 ' 



x( 



^- x ( v + i Q)^ x ( v -iQ) 



x(v + ij) 
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They satisfy the following relations 

Kf(u, v)= l - (K yQ {u, v) =f K Q {u - v) 

i d % % 

K vq(u,v) = — — log S yQ (u,v) = K(u,v + -Q) - K(u,v Q) 

2m du g g 

x(u) - x(v + iQ) -7-t - x(v - i®-) 
S yQ (u,v) = S^(u,v)Sf(u,v) = .q — p 

x[u) — x[v — t—j -rr — X [V + I— , 

S(u,v + IQ)x(v-tZ) 
~ S (u,v-±Q)x(v+i2) " 

It is worth mentioning that S v q(±2,v) = 1. 

Next, we define the following S-matrices and kernels 

x(u — i—) — x(v + i—) x(u — i—) — x(v — i—) x(u + i— 
Sg M {u,v)- K " ;; "' •" 9 



x{u + £— ) — x(v + iy) x(u + i— ) — x(f — Zy) x(u — ij 

M-l „, „, i 



and 



V g ! \ g I \ g 1 \ > g ! \ g ) 



3 



A. u - v + ±( M -Q + 2j) 



(8.9) 



nu-v- % ~{Q-M + 2j) 
^ 4, (8-10) 
j= i u-v + ±(Q-M + 2j)' K } 



~ x(u- i 3 -) - x(v - i—) x(u + i 9 -) - x(v - i K ) x(v + i—) 

\ g ) \ g ) \ I g ) \ g ) \ g ) 

^u-v- i(M-Q + 2j) 

x rr ii — z — iL t (s.ii) 



Then, the s\{2) S-matrix in the uniform light-cone gauge |59| with the gauge 
parameter a = can be written in the form 

= S^iu-vr'ZQuM- 2 , (8.12) 

where Tfi M is the improved dressing factor JBE ]. The corresponding $1(2) and dressing 
kernels are defined as usual 

K S(?)( U , v ) = ^--^ lo S S %™)( u , w ) ' Kqm( u i v ) = ^-^ log^QAfK • (8-13) 
The analytically continued st(2) S-matrix is given by 



S ^ V) S lM (u-v)Z UM (u,vy 
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where the improved dressing factor is given by 

- log S UM (u, v) = $(y+, y+) - y 2 ) - Qfc, y+) + $(yf , ) 

1 



+ g O^W) + ,Vi) - ^(2/^,2/1 ) - HV2,Vlj) (g 14) 

i (yi-yt)(yi-±)(yt-±) 

■ — log- 



2i "° {yi-yi)(yl-j,) 2 

Here = x s (m ± |) are parameters of a fundamental particle in string theory, and 
vt = x ( v i g^O are parameters of an M-particle bound state in the mirror theory, 
see 16131 for details. 



Next, we introduce the following kernel and S-matrix 



K ( u > v ) = 7T~ T" lo S Sma(u, V) = ^ " — — , S mfl (u, u) = ^\ . 

2m du 2n V4 - u 2 u - v x(u) - ^-^y 

With the help of this kernel we define 10 

K(u,v) = K(u,v)[9(-v - 2) + 9(v - 2)] , (8.15) 

K Q {u,v) = [K{u+-Q,v) + K{u--Q,v)][6{-v-2) + 6{v-2j] , (8.16) 
9 9 

where 6{u) is the standard unit step function. Obviously, both K and Kq vanish for 
v being in the interval (—2, 2) and are equal to (twice) the jump discontinuity of the 
kernels K and Kq v across the real semi-lines \v\ > 2. 
We also use 



1 d , „ , , 1 V 1 ^ v „ , s sc. (it) — xJv) 



K ss (u,v) = — — log S ss (u,v) = — . = 7 r, S ss (u,v 



2ni du 2m L T(u — v)' x s (u) 



l ' 



and define Sq(u,v) as 

^(■u — -Q, v) I d - 

Sn{u,v) = ; f -, Kn(u,v) = — log So (u, V ) , (8.17) 

QK ' ; S ss {u+ - g Q,v) QK J 2m du & QK ' v ; 

to ensure that Sq(— oo,v) = 1. 
The quantity £ is defined as 

£(u) = hg-) f = 21ogx s (w) ^ for u 6 (-oo, -2) U (2, oo) . (8.18) 

x[u + 10) 



10 



The definitions of the kernels K and Kq differ by the sign from the ones used in 22 
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The TBA equations for Yg-particles involve the kernel 



k ° = ^i~L hgtQ = - K ^^ + ^ (8 - 19) 



where 

oo „2 



I Q = Y^K 2n+Q {u ) v) = Kf +2 \u-v) + 2 dtK [ ® +2] (u-t)K(t,v), (8.20) 
n=i J - 2 

^ r (M) = ^ log r[fl + >] = ^ + S ^ 2n+Q - 2(u) " • (8 - 21) 

The kernel ( |8.19| ) is related to the dressing kernel Kq M as follows 

K%(u, v) = K^(u, v+ Z --iO) + K% 1 (u,v-^ + zO) - K* 2 (u, v) . (8.22) 
The analytically continued kernel Kf^ is given by 

Kl (u, v) = -K Uy iI -K„(u--,v). (8.23) 
Finally, integrating ( |S.23| ) over the first argument, one gets 

]ogt u (u,v) = - log S Uy * I - log S aa (u- -,v). (8.24) 

Note that the first term here is real, but the last one is not. 

8.2 Solution of the Bethe-Yang equation for the Konishi state 
Perturbative solution 

The equation (|3.2fl can be solved in perturbative theory, and one gets up to g 16 



_ 2ix V3g 2 , Qy^ff 4 3^ 6 g 8 v / 3(960C(3) + 960C(5) + 671) 

P Y " ~T~ + ~2 8~ 9 (C(3) + 1} + 1024 (8 " 25j 



3 10 



r- , 141C(3) _ 309C(5) _ 315C(7) _ 3807 \ 
V 64 128 128 2560/ 

/ 2799C(3) 9C(3) 2 1527C(5) 3339C(7) 441((9) 7929 \ 
V " 1 512 64 256 512 64 2048 / 

30015C(3) 81C(3) 2 7929C(5) 45C(3)C(5) 17127C(7) 



1197C(9) 10395C(H) 303837\ 
64 512 28672 J 

16 ^/ 340785C(3) 2349C(3) 2 350505C(5) 891C(3)C(5) 225C(5) 2 
9 V 8192 1024 8192 256 256 

183519C(7) 945C(3)C(7) 100863C(9) 230175C(H) 127413C(13) 15543873\ 
+ 4096 + 512 + 2048 + 4096 + 2048 + 524288 J ' 
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Approximately one gets 

p = 2.0944 - 0.433013/ + 0.487139/ - 1.43028/ + 4.77062/ (8.26) 

-15.7964/° + 52.5014/ 2 - 176.638/ 4 + 602.45/ 6 . 

The corresponding expansion of the w-variable is given by 



1 f l + 2/ - b 4-+9 6 (- A + ^P) t^.9 8 (461 + 144C(3) + 240C(5)) (8.27) 



u = 



V3g 

gl0 /1133 + 63C(3) + 189C(5) + 315((7)\ 



4 * \A 4 J 128- 

) 



.12 



128 32 64 64 

( 23835 1167C(3) 9((3) 2 729((5) 2079C(7) 441((9)\ 



9 V 1024 + 256 + 64 + 128 + 256 + 32 ) 
14 / 64731 3429C(3) 9((3) 2 897((5) 45 8559C(7) 
+ 5 v l024 + 256 + ~128~ + ~~ 64~ + 64 C(djC(5) + 512 

189C(9) 10395<(11)\ 
+ 8 + 256 J 
16 / 1441077 176445C(3) 405C(3) 2 169371C(5) 477 
9 \ 8192 4096 512 4096 512^ M ' 

225((5) 2 87417C(7) 945 51975C(9) 147015C(11) 127413C(13) \ 

256 2048 512^ ' 1021 2048 1021 J 

Approximately one gets 

57735 

u = - + 1.15470 - 0.721688/ + 1.53087/ - 3.98263/ (8.28) 

9 

+ 11.1101/ - 32.8297/ 1 + 101.602/ 3 - 325.587/ 5 . 
The dimension of the Konishi operator is 

A = 2 + 3 9 2 - 3/ + ^ + / f-^ - « V , 10 f ^ + + ««5)\ 

y y 4 y V 64 8 / V 256 32 16 y 1 ' 



+ 9 12 



y v 2048 256 128 512 256 64 y 

16 / 7964283 _ 278505C(3) _ 621C(3) 2 _ 58491C(5) _ 135 

9 V 16384 2048 512 512 128 

198207C(7) 20979C(9) 31185C(11) 
2048 256 512 

/22613385 3600585C(3) 1539C(3) 2 1520127C(5) 7101C(3)C(5) 



67287 27((3) 1377<C(5) 945<C(7)\ 
1024 2 128 128 J 

Z359655 10899C(3) 27C(3) 2 18117C(5) 7371C(7) 1323<C(9)\ 



+ 9 16384 + 8192 + 256 + 4096 1024 

675C(5) 2 2605095C(7) 2835<C(3)C(7) 573237C(9) 1002375C(H) 382239C(13)\ 
+ 512 + 8192 + 1024 + 2048 + 4096 + 2048 J ' 

Approximately one gets 

A K onishi = 2 + 3/ - 3/ + 5.25/ - 12.3679/ + 33.8743/° - 100.537/ 2 (8.30) 
+ 313.532/ 4 - 1011. 73/ 6 + 3348.ll/ 8 . 
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Figure 11: Numerical solution of the BYE for the Konishi operator. The momentum was 
computed for < g < 10 with the step Ag = jq, and for 10 < g < 100 with Ag = 1. On 
the right picture the asymptotic dimension of the Konishi operator is plotted. It approaches 
2^2ng as expected from [ p6| . 



We have also solved the BY equation numerically for small values of g, and its 
numerical solution perfectly agrees with the analytic one. The perturbative solution 
works very well at least up to g — |. For g = \ the difference between the analytic 
solution and numerical one is ~ 5 x 10~ 10 . 

Numerical solution for i < g < 100 

The BY equation can be solved numerically up to very large values of g, and one 
gets the plot on Figure 7. for the momentum p as a function of g. For large values 
of g ~ 100 the momentum is approximated by 

/2ir 1 , 

Pafs = \ , 8.31 

V 9 9 

with a good precision as expected from |56|, |32|. For g = 100 the difference between 



the numerical solution and the AFS formula is equal to —0.0016902. If one uses the 



following asymptotic expression from [58 



2tt 1 0.931115 + 0.199472 \og(g) 




Prs = \ + ^ , (8.32) 

V 9 9 9 3/2 

that is one includes the next subleading terms then the agreement is much better and 
for g = 100 the difference between the numerical solution and the Rej-Spill formula 
is equal to 0.0001587. To match the coefficients in p-^s one has to solve the BY 
equation for larger values of g. 

The corresponding plots of the u- variable are shown on Figure 3, and were dis- 
cussed in section 3. 

8.3 Transfer matrices and asymptotic Y-functions 

Here we remind the construction of the asymptotic Y-functions in terms of transfer- 
matrices corresponding to various representations of the centrally extended su(2|2) 



superalgebra. Consider K 1 physical particles (excited states) of string theory char- 
acterized by the w*-plane rapidities u±, . . . ,un, N = K 1 or, equivalently, by the 
Pi, . . . ,pn physical momenta. Each of these particles transforms in the fundamental 
representation of su(2|2). Consider also a single auxiliary particle with rapidity v 
corresponding to a bound state (atypical) representation of su(2|2) with the bound 
state number a. Scattering the bound state representation through the chain of iV 
physical particles gives rise to the following monodromy matrix 

N 



T(v\u) = Y\_$oi(v,Ui) 



i=i 



Here S>oi(v, Ui) is the S-matrix describing the scattering of the auxiliary particle in the 
bound state representation with a physical particle with rapidity itj. The transfer- 
matrix T a (v\u) corresponding to this scattering process is defined as the trace of 
T(v \u) over the auxiliary space of the a-particle bound state representation T a {v\u) = 
tr T(v\u). We are mostly interested in the situation where the auxiliary particle is 
in the mirror theory. An eigenvalue of this transfer matrix for an anti-symmetric 
bound state representation of the mirror particle is given by the formula found in 



T a ( v \u)=n 



Di—x+ V x 



i=l 




K 



V- TT yj-x- 
L i 11 Vi-x+ 



i=l 



a-1 K n 

\~* TT n. •■ / •'• 

L. 4 11 Vi-x 
k=0 i=l 



i=i 

K 1 

n 



i - 



2ik 
9 



v- Ul + Ua-l) 



i=l 



X 



x± 9 i i i y_i 

^1^11 tt + i-x+- 4 



n 



8=1 



K III 



i 1 . i(2k-\ 

Wi—X + V + — — ■ 

1 x+ g 



+ 1 . i(2fc+3) 



8=1 



2i^+l) H m _ x +_J_ + i(2k+l) 
i=l 



(8.33) 



Eigenvalues are parametrized by solutions of the auxiliary Bethe equations: 




TT _ yK Vk g 

i=i 1 UK Vk g 

jt «'*-«'< + 7 



i=l,i^fe 



1 g 



(8.34) 
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In the formulae above 



v — x 



1 % 

x H 1 — a , 

x- g 



x^ g 

and v takes values in the mirror theory v-plane, so x^ = x(v ± -a) where x(v) is the 
mirror theory x-function. As was mentioned above, Uj take values in string theory 
it-plane, and therefore xf = x s (uj ± ~) where x s (u) is the string theory x-function. 
Finally, the quantities A± are 



\±(v,Ui,k) 



[x, x 



1)(XH 



X 



+ 



x + (x,- + xj 



(x i — x+)(x+x+ — 1) g (x- — x + )(x+x 
ix + (x7 — xf) 



1) 



(x, — X+)(x + Xj 




i(2k - a) 
9 



(8.35) 



For the sl(2)-sector X 11 = = K and the expression above simplifies to 

K 1 



T a (v I u) 



(x — x i )(1 — X X+ ) X 4 

1 + -srK 1 - 



a-1 K L 



- 2 En 




X+ — X 



X^X,- X" 



1 



(8.36) 



2ik 
9 



Ui + ±{a 



1; 



a-l if 1 



2^ 2^|| A m (t7,Mi, fc) 
m=± fe=l i=l 



fc=0 i=l 

For a two-particle physical state, the last formula appears to coincides up to a gauge 
transformation with p8" 



T a (t> |u) = — 
x~ 



X 



x„- 



+ a 



i=l 
2 

n 

i=i 



X" 1 " — Xj 

x + — x7 



z7 



i=l 



nr~v~ Of rp ! 



n 



X 



x, x,- — 



1=1 1 1 



1 

:r + 



(8.37) 



which is nothing else but an eigenvalue of the transfer matrix evaluated on the 
fermionic vacuum and continued to the mirror region for the auxiliary variable. In 
numerical computations formula (|8.37 ) works much faster than ( |S.36| ) and, therefore, 
we use it for constructing the asymptotic Y-functions. 

The transfer matrices T ayS are used to introduce Y-functions Y a)S which solve 



Y-system equations in the following standard way |33], |34|, |27 

Tq, s -L.±T a s— 1 



T a +i,sT a -i, 



.38) 



The transfer matrices T a s solve the so-called Hirota equations |j0] and they can be 
computed via T a l = T a {y \ u) through the Bazhanov-Reshetikhin formula which 
we present here for the asymptotic solution 



T a , s = det 1 <ij< s T a+i _ j [v + -(s 



J) I u 



(8.39) 
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If here a + i—j < then the corresponding element T a+i „j yu + ~(s + 1 — 

regarded as zero. Also, T +^(s + l— i— j) \ uj = 1. In other words, asymptotic T Qj 
and T ls satisfy the boundary conditions: T 0jS = 1 = T a , T a <o,s = and T , a s<0 = 0. 

Here for reader's convenience we recall the relationship of the Y-functions intro- 



j) | u) is 
i 



duced in [22] and those of \27\. For the auxiliary Y-functions we have 



y(l) = 


1 


^2,1-?0,1 


y(2) 


1 


T2, -1^0,-1 




^1,2^1,0 ' 




^1,0^1,-2 


yW _ 


^2,2 — 


^2, 3^2,1 


y(2) 


— ^2,-2 — 


^2,-1^2,-3 


^3,2^1,2 


^3,-2^1,-2 


= 

Q|mt> 


1 


7g+2,i^g,i 


y(2) 


1 


Tq+2,-iTq-X 




^Q+l,2^Q+l,0 




Tq+1,oTq + 1-2 


^ = 


Yi,q+i — 


Tl,Q+2T\ t Q 






T\-qTi-Q-2 


^2,(5+1^0^+1 





140) 



The asymptotic functions Yq{v) in the sl(2)-sector are given by ( |3.8| ). Note that in 
this expression the transfer matrix naturally comes with a prefactor 



K 1 K 1 



Sq(v \u) = n y/SSfoiv, Ui) = ft , 1 I (8-41) 



, - ! »=i \/S e ^{Ui,v) 



This prefactor can be split as follows 



K 1 _ if 1 



Sq(^ I n) = n , 1 = II \l h Q^ v ) » ( 8 - 42 ) 

i=1 y/S s i^(ui,v)hQ{ui,v) k =i 

where }iq(u,v) is introduced in ( |S.97| ). This quantity satisfies the equation 

h Q (u, v + -)h Q (u, v - -J = /iQ+i(w, v)h Q -i(u, v) 

and is a pure phase when u and w are in the string and mirror regions, respectively. 
The splitting (|8.42|) can be used to introduce the normalized transfer-matrix Tq i (v) 



K 



Tq,i(v \u) = ]JJ h Q (u k , v) T Q (v\u 



k=i 



which renders the corresponding Yq = e J£q Tq i(v\u) real. The functions Yq repre- 
sent a simple and useful tool for checking the TBA equations. 
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8.4 Hybrid equations for Kg-functions 

In this appendix we derive the hybrid TBA equations for Ig-functions. We discuss 
only the ground state TBA equations because the excited states equations can be 
obtained by using the contour deformation trick. 

The canonical TBA equation for Q-particles reads 

log Y Q = -L£ q + log (1 + Y Q ,) * K%% (8.43) 
+ 21og (l + — L-) *K%? + 21og (l - i-) iK* + 21og (l - -L) IKf . 

To derive the hybrid equation we need to compute the first sum on the second line 
of this equation. It is done by using the simplified equations for Y M \ VW 

\ogY M \ vw = log(l + *Af-i|«tB)(l + Ym+i\vw) * 8 (8.44) 

i - y_ 

- log(l + Y M+l ) *s + 5 M i log — * s. 

1 — Y + 

Let us assume that we have some kernels Km which satisfy the following identities 

Km-s* (K m+ i + K M -i) = 5K M (M > 2) , K x - s * K 2 = 5K, , (8.45) 

where SKm are known kernels. Then, applying the kernel Km to both sides of (|8.44j) , 
and taking the sum over M from 1 to oo, we obtain the formula 



( 1 + 77 ) * K A/ = ^ l0g i 1 + Fm I^) * 5Ka/ 

M=\ ^ M=1 



146) 



i-y_ 



+ X log(l + Ym+i) * s * K M - log — — * K 



M=l 



Now choosing the kernel Km to be Kf£® and using the formula (6.31) from W2 



K%2 ~ ' * [KZf' Q + K^ Q ) = <W,Q * (M > 2) , 

K% x - s * ^ = s * AT yQ + 5 2iQ s , (8.47) 

and the identity (|8.46|) , we get 

V log f 1 + j * K^Q = log (l + Y llvw ) *s*K yQ + 8 M +i,Q log (l + Y M \vw) * s 

oo 1 - y 

+ X log(l + F M+ i) * s * - log I — ^ * s * ^ . (8.48) 



M=l 



Finally, substituting this formula into ( B.43Q , we get the hybrid ground state TBA 
equation. 
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8.5 Analytic continuation of Yi(v) 

To derive the exact Bethe equations (|4.17|) one has to analytically continue Y\{z) in 
eg . (|4 . 1 3|) to the point z^. On the f-plane it means that we go from the real f-line 



down below the line with Im(t>) = — ^ without crossing any cut, then turn back, 
cross the cut with Im(t>) = — and |Re(t>)| > 2, and go back to the real v-line. As a 
result we should make the following replacements x{v — |) — > x s (v — ^) = x(v — p, 
x{y + i ) — > x s (v + -) = l/x(v — 1 ) in the kernels appearing in (|377|) . 

The analytic continuation depends on the analytic properties of the kernels and 
Y-functions. Let us consider the terms in eq.( |4.13| ) order by order 

Terms with Kj_(t>)-functions 

For the analytic continuation of the last two terms in ( [4.13 ) it is convenient to use 
the kernels K± given by 



1 1 d i i x ( u ) ~ x ( v + g) 

K^L (u,v) = - — :— log 5*^ (u, v ) , S* y (u,v) 



x(u) — x(v — ~) ^ x(v + ~) 



x(v — -, 

v 9' 



\ d -^r — x(v — - 

K\ (u,v) = — — logS^ (u,v) , Sl(u,v) = ^- 



<j' 



x(u) V 9 7 \ 



_J 9j_ 

x(v — -) 



.49) 



(8.50) 



It is clear from these equations that K^(u,v) remains regular in the analytic con- 
tinuation. K v __ (u } v) on the other hand has a pole at v = u — - and behaves as 



11 i 

K y _}{u,v) = T + regular at v ~ u . (8.51) 

2ni u — v — - ' q 

9 a 

Thus, one needs to analyze the analytic continuation of a function defined by the 
following integral for real v 

Hv) = ^ f du f(u) l — . (8.52) 

2,1X1 J_2 u — v — - 

9 



The consideration is the same as the one for the dressing phase in [66], and after v 
crosses the interval [—2 — ~, 2 — ~] one gets the following expression for F{v) 

F(v) = f{v + -) + / d M /(«) - , Im(i;) < -- . (8.53) 

g 2m J_ 2 u-v-j g 

Then we go back to the real f-line but we do not cross the interval [—2 — -, 2 — -1, 
and therefore ( |8.53| ) remains the same. However, we should also analytically continue 
f(v + -) back to real values of v. Thus we conclude that the analytic continuation of 



log ( 1 - — 1 *K_ + log ( 1 — — J *K? (8.54) 
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is given by 



log ( i - y-^ry j + log (i - ^) +iog (i - ± ) tKf* , (8.55) 



where Y*-(v) is the Y-(v) -function analytically continued to the upper half-plane 
through the cut \v\ > 2. 

As was discussed in section 4, coincides with Y + , and, therefore, to find 
Y*-(v + |) one can just use the analytic continuation of the TBA equation QS.69Q 
for Y + {y) to Y + (v + 1 ). Since the S-matrix Sj v (uk,v) has a pole at v = Uk + ^, see 



appendix |Q| , one concludes that Y*_(uk + -) = oo, and, therefore, the contribution 
of the first term in (|8.55|) vanishes in the exact Bethe equation Yx„{uk) = —1. 

Nothing dangerous happens with the term log * s * K^ x because there is 
no singularity in the analytic continuation process of K^ wx until we get back to the 
real w-line. Then, we get for the analytically continued K^ x 

K lwx( u > v ) = 7T~ ~T~ 1°S S lwx(u, v) , (8.56) 



2ni du 

^vwx{ U l v ) 



i 

11. („, „,\ — " v ~ 9' ' 9' ^"~a 



x(u-±)-x a (v + ±) W < r 



x ( u+ i)- Xs ( v -i) -^- Xs ( v + i) 



and it shows a pole at v — u 

1 i 

Kill (u, v) = h regular at v ~ u . (8.57) 

2%i u — v 

Since we integrate over the line which is above the real line, the pole is not dangerous 
if the function we integrate with the kernel is regular as it is the case for the case 
under consideration. 

Terms with Ym ^-functions 

The analytic continuation of the term log (l + Yi\ vw ) * s*K y i is given by 

\og(l + Y 1{vw )*(s*K yl + s) (8.58) 
because K y i{u, v) has a pole at v = u — K 

Terms with YQ-functions 

The kernel K^Ls ^ s gi ven by 

= -K Q1 (u-v)-2K% 1 (u,v) (8.59) 
= -Kq-i(u -v)- K q+1 (u -v)- 2^ (u, v) , 
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Since Kq^u, v) is a holomorphic function if u is in the mirror region and v is in 
the string one, only Kq\ can cause any problem with the analytic continuation. 
Moreover, we see immediately that only the Q = 2 case should be treated with care. 
It is easy to show however that since the integral over u is taken from — oo to oo 
the analytic continuation does not give any extra term because v crosses the line 
Im(w) — ~g twice. This is the difference of this case from the Y± one. The term 
with K^~ l,l {v) is harmless too because s is regular. 

The term log S * K^^uJ ,v) 

This term is similar to log ( 1 + ^J— ) * KVl considered above, and its analytic 
continuation is given by 

N . x - L 

log Res(S) * K^ x (uj, u k ) - 10 § («i -«*--) ~~~ ~r ■ ( 8 - 6 °) 

,=1 3 X j X k 



Summing up all the contributions, one gets the exact Bethe equations ( |4.17 ) and 
( 3.72 ). We have shown that the r.h.s. of (|4.17 ) is purely imaginary, and, therefore, 
the exact Bethe equations can be also written in the form 



A' 



ix(2n k + 1) = Lp k + i ^ lc, g ^(2) Uk ) ( 8 - 61 ) 



2k x J 



N N 



+ 2 V log Res(S) ★ lmKH* x (uj } u k ) - 2 V Im log ( Uj - u k ) 

~t ~t 9 x j ~ x k 

- 2 log (1 + Y Q ) * (lmK% u - s * 1mK%£> u ) - 2% log (l + Y 1{vw ) * (§ + s*K vU ) 

- 2 log ^-y- * s * ImKUl - i log (1 - — ) (1 - i-) *K vU . 

Analytic continuation of the canonical TBA equation for Y\ 

Let us also consider the analytic continuation of the canonical TBA equation for Y\. 
Then, the kernel K^ x is given by 

KHM = ±-*-]o S S™(u,v) , (8.62) 
qMi i , ^-jM)-,(, + i)^-^i) 

(It, V 



VWXK ' ' + |Af )-«;(«- |) - ^(t, + ±) ' 

and nothing dangerous happens for M > 1. So, we just have to consider the M = 1 
case. Since Yi\ vw (uj) = 0, and we should take care of the resulting log-singularity 
before the analytic continuation. 
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Introducing Z-functions as in flS.73Q , one gets for the term in the canonical TBA 
equation ( p.68|) 



2 log 1 + 



V 



M\vw 



M 2 



N 



-11 



M 2 



- * KZ' 1 ' 1 - 2 log J] (u - Uj ) * iC, 



3=1 



N 



+ 2 log Z M \vw * K, 
2i, x 



Ail 



- 2 log 2 - 2 ^ log ( Uj - v - -) i_Y + 2 log Z M| ™ * 
J= i 51 x i z+ 

where we used log tMzi * ^muz 1 ' 1 = 1°S 2, and the following formula 



(8.63) 



dt log(t + «0 - iC*(*, u) = log - v - -) ^ — 2f 



y T 



The analytic continuation in v then gives 



Ml v w 



N 



2 log ( 1 + ^— ) * lfj£r(«) = 21og2 - 2^>g («, - « - -) 



2i x - 



J a;- 



g Xj — ar 



+ 2 log Z M]vw * p .„. + log Z l \ vw {v) . 



The equation seems to coincide with the one derived in |3(J after one performs the 
parity transformation Uj — > —Uj, and sets g = 2. 



8.6 Canonical TBA equations 
Canonical TBA equations: g < g£) 



The canonical ground state TBA equations |22| are derived from the string hypothesis 
for the mirror model [21] by following a textbook route, see e.g. [j53 |, and can be 
written in the form 

• Q-particles 



Vn = loeFo + LSo - log (1 + Y a >) * K 



•1(2) 



21og 1 



jy-M'Q 



-21og(l-^) */^-21og(l-i_) = 0, (8. 



64) 



^-particles 



V± = lo g y± + log (1 + Kg) * Kg" - log 



1 + 



1 + 



1 



* K M = . 



(8.65) 
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M | vw- strings 



V M \vw = log Y M \ VW + log (1 + Y Q ,) * if£ M 

- log (l + ■—— ^ * ^-m'm - log ^ * Km = . (8.66) 

V *M'\vw J YZ 



• M|w-strings 

( 1 \ 1-^- 
Viif |w = log Ym\w - log 1 + — * K M , M - log r *K M = . (8.67) 

Applying the contour deformation trick to the canonical TBA equations, one 
gets the following set of integral equations for Konishi-like states and g < g& 

• Q-particles : Vq — - ^ log S£* ( S (Uj, v) , (8.68) 

• particles : V± = / ] log S± y (iij, v) , (8.69) 

• M|tw-strings : V M \vw = ^ log S^ M (u h v) , (8.70) 

• M|w-strings : V M \ W = . (8.71) 

Here summation over repeated indices is assumed. The sums in the formulae run 
over the set of particles, all Y-functions depend on the real u variable (or z) of 
the mirror region. We recall that S^S and S],* V M denote the S-matrices with the 
first and second arguments in the string and mirror regions, respectively, and both 
arguments of the kernels in these formulae are in the mirror region. 

Then, the integrals are taken over the interval [—2, 2] for convolutions involving 
Y±, and over the horizontal line that is a little bit above the real u line (or the interval 
Re(^) G (— *y), lm(z) = on the 2-torus) for all other convolutions. 

The reason why one should choose the integration contour to run a little bit above 
the real line of the mirror w-plane is that Yx\ vw has zeros at u = Uk, and, therefore, 

the terms log ( 1 + -r^- — ) * K with any kernel K should be treated carefully. The 
prescription above for the integration contour guarantees the reality of all Y-functions 
as we show in appendix |3.7| . The log S-term in the equation for vw-strings is in fact 
necessary to cancel the corresponding singularity on the l.h.s. of this equation. 

One can show that the imaginary zeros of Yk\ vw and 1 + Yu vw do not contribute 



to the canonical equations for Konishi-like states at weak coupling. In appendix fTB 
we also show that the simplified equations can be derived from the canonical ones 
following the same route as in ^2), p3| . 

The canonical TBA equations are rather complicated and involve infinite sums 
that makes the high-precision numerical tests very time-consuming. We have checked 
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numerically that for Konishi-like states they are solved at the large L limit by the 
asymptotic Y-f unctions if L = J + 2. 

Let us stress again that the TBA equations above are valid only up to the first 
critical value of g where the function Yi\ vw {u) has real zeros only for u = u k , and 
^MiDw-functions with M > 1 do not have other zeros in the physical strip. 

Exact Bethe equations: g < gil' 

To derive the exact Bethe equations we take the logarithm of eq. ( |3.7| ), and analyt- 
ically continue the variable z of Y\{z) in eq.( |B.68|) to the point z* k . The analytic 
continuation is similar to the one in section |], and its detailed consideration can be 
found in appendix |3.5| . As shown there, the resulting exact Bethe equations for a 
string theory state from the sl(2) sector can be cast into the following integral form 

ni(2n k + 1) = logY u {u k ) = iLp k - X) lo g 5 i(2)( u i> u k) + log (1 + Y Q ) * K® 1 



.1(2) 



2i^ x j 



+ 2 log 2 - 2 V log ( Uj -u k --) — x -\ + 2 log Z M[VW * Kf± 

g Xj - x+ 

+ 2 log (l - i-) *K1 U + 2 log (l - i-) *K y + u . (8.72) 

Here the integration contours run a little bit above the Bethe roots Uj, p k = iSq(z^) = 

—i log —, jt is the momentum of the fc-th particle, and the second argument in all 

the kernels in this equation is equal to u k but the first argument we integrate with 
respect to is the original one in the mirror region. The Z-functions are defined in the 
same way as in [[SIJ 

1 4 1 1 (M+l) 2 

1 + iv " Zllvw 3 nf =1 (u - Uj) ' 1 + y^z' Zm1vw m(m+ 2 )' (8 - 73) 

they are positive for real w, and Zm\vw, M = 2, 3, . . . asymptote to 1 at u — > oo. 
Taking into account that the BY equations for the sl(2) sector have the form 

N 

m(2n k + 1) = i J p k - ^ 1°S S li\l) ( u ji u k) , 

3=1 

and that Yq is exponentially small at large J, we conclude that if the analytic con- 
tinuation has been done correctly then up to an integer multiple of 2ni the following 
identities between the asymptotic Y-functions should hold 

A 2i X J-^ 

TZ k = 2i Pk + 21og2 - 2^1og ( Uj -u k --) _ x \ +2\ogZ M \ vw *KZ 1 * 

j=l 9 X j X k 

+ 21og(l - ^)*K y _ u + 21og(l- ^r)*Kf* =0. (8.74) 
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For N = 2 and U\ = —u 2 one gets one equation, and by using the expressions for the 
Y-functions from appendix one can check numerically that it does hold for any 
real value of U\ such that only Y\\ vw has two zeros for real u. 

We believe that the canonical TBA equations (|8.68| - |8.71 ) and eqs. (|3.72D are 



equivalent to the ones proposed in |30] for states having the vanishing total mo- 



mentum but a detailed comparison is hard to perform due to different notations and 
conventions. We see however that they are definitely different for physical states 
which satisfy the level-matching condition but do not have the vanishing total mo- 
mentum. 

Canonical TBA equations: g^r < g < g$ 

Here is the set of canonical TBA equations for Konishi-like states and g$ < g < g~cr 

2 

. Q-particles : V Q = - £ log {u 3 h v) - 2 ^togS^-*) (8.75) 

* i=i 

+21og^(r 1 , U ) + 21og^ Q (r 1 , U ), 

N 2 

• y— particles : V± = log S± y (Uj, v) — log S\{rJ — v) + log S^n — v) , 

3=1 3=1 

N 2 

• M|fw-strings : V M \ VW = ^log S^* M (uj , v) - ^\og S 1M (rJ ,v) + \ogS 2M (r 1 ,v) , 

3=1 3=1 

• M\w— strings : Vm\w — , 

where 2logS y Q(rx, v) appears due to the imaginary zero r\ of Y±, and we take into 
account that the S-matrix S v q is normalized as S y Q(±2,v) = 1. 

Exact Bethe equations: g$ < g < g^ 

The exact Bethe equations are obtained by analytically continuing logYi in QS.76] ) 
following the same route as for the small g case, and they take the following form 



N 

QU 

1(2) 

i=i 



ni(2n k + 1) = log Y u (u k ) = iLp k -^2 lo S S li{l) u k) + log (1 + Y o) * K > 



- 2 log S] l w ' x {r j , u k ) + 2 log SZlin, u k ) + 2 log S yU (n,u k ) 

3=1 

N . X T L 

- 2 log («i -u k --) = _ x \ + 2 log 2 + 2 log Z M]VW * KZl* 

j=l 9 X j X k 

+ 2 log (l - ^) ZKl 1 * + 2 log (l - -L) icKf* . (8.76) 
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We conclude again that the consistency with the BY equations requires the fulfillment 
of the identities between the asymptotic Y-functions similar to (|8.74|) that we have 
checked numerically for the Konishi-like states. 



Canonical TBA equations: g$ < g < gi? 



The canonical TBA equations for Konishi-like states and gil' < g < gi? take the 
form 

N 2 

• Q-particles : V Q = - ]T log S^u,, v) - 2 ]T log S^rj ,v) , (8.77) 

3=1 3=1 
N 2 

• y— particles : V± = log S± v (v,j, v) — log S\(rJ — v) , 

3=1 3=1 

N 2 

• M|mu-strings : V M \ VW = log S^ M (uj, v) - ^\og S 1M (rJ ,v) , 

3=1 3=1 

• M\w— strings : Vm\ w = , 

The equations for Yg-particles differ from (|8.76| ) by the absence of the term 2 log S%® x (ri, v) 
The p.v. prescription however would give instead additional terms due to the real 
zeros Tj of Y2\ vw and Y±. 

Exact Bethe equations: g& < g < gi? 



Analytically continuing log Y\ in (|8.77|) , one gets the exact Bethe equations 

N 

m(2n k + 1) = log Y u (u k ) = iLp k -^ lo S S li(l) ( M i> u k) + lo S i 1 + Y o) * K t 



N 

Ql* 

1(2) 
3=1 



2L x 3~- 



- 2 J] log SltirJ, u k ) + 2\og2-2j2 log («j -«*--) 

3=1 3=1 9 X i 

+ 2 log Z M \ VW * K™1* + 2 log (l - i-) * A! 1 * + 2 log (l - ^) * K y + U . (8.78) 

We recall that the integration contours should run a little bit above the Bethe roots 
Uj, and below the dynamical roots rj, and the consistency with the BY equations 
leads to identities of the form ( |8.74| ) that we have checked numerically. 

Canonical TBA equations: g& < g < g^f 1 

The necessary modification of the integration contour follows the one for the first 
critical region, and the integration contour runs a little bit above the Bethe roots Uj, 
below the zeros rj , k — 1, . . . , m, and between the zeros rj m+1 \ Im(r[ m+1 ^) < 0. All 

the roots — - are between the integration contour and the real line of the mirror 

3 g to 

region. The contour for Y± functions lies above the zeros r . 
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Here are the canonical TBA equations for Konishi-like states and gffl < g < gffl 

N m-1 2 

. g-particles : V Q = - £ log S^fav) - 2 £ £ log S^rf v) 

3=1 k=l j=l 

qm~l,Q( r (.m+l) \ q m+1 Q ( (m+1) \ 

,01 J vwx VI i u °vwx VI i U ) /a vnN 

+2 log — r — -r t , (8.79) 

*~>vwx V \ i u ) 

N m-1 2 

• y— particles : V± = log S± v (v,j, v) — log S k (rj k — v) 

j=l k=l j=l 

S m -i(r { r +1) ~ v)S m+1 (r^ +1) - v) 



+ log 



„(m+l)- \ ' 



N m-1 2 

• MK-strings : V M]vw = ^ogS 1 x ; M (u„v) - ^ ^ log S kM {rf 

i=i fe=i i=i 

c f r ,( m + 1 ) „,\ c /V ( m+1 ) ,,1 

, Jm-lM^l i^JJm+l.Ml^i ,V) 

"m,M V 1 ' y / 

• M\w— strings : Vjw| w = , 

Exact Bethe equations: gffl < g < gffl 

The exact Bethe equations take the following form 

N 



t(2n k + 1) = log Y u (u k ) =tL Pk -J2 log Sj,*(J) u *) + 1°S ( x + *q) * K $w 

3=1 

m 2 ™-l,l,/„Kl) ,, \ qm+1,1, { Am+l) \ 

- 2 ^ ^ log ^;(rf^- u k) + 2 log g - ( \ m1 :;jifir 



fc=i j=i 



Svwx( r l i U k) 



2L 



+ 2 log 2 - 2 V log («j - Wit ) — -qr + 2 log Z M | W * iC 

j=i y j fc 



Ml. 



+ 2 log (l - i-) + 2 log (l - i-) **T?* , (8.80) 

where xj = x(rj p, and the integration contours run just above the Bethe roots, 

and below the zeros A , k = 1, . . . , m. The functions Z M \ VW are defined as in (|8.73|) . 

The consistency with the BY equations again implies that the sum of the terms 
on the last three lines of (|8.80|) is equal to —2ip k that we have checked numerically 
for the Konishi-like states. 

Canonical TBA equations: g~cr < g < gd? +l ' 

(k) 

In this case all the zeros , k = 1, . . . ,m + 1 are real, the integration contour 
runs below all of them but the Bethe roots Uj, and the canonical TBA equations for 
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Konishi-like states and gi™ < g < take the following form 

N m 2 

. Q-particles : V Q = - £ log S^fa, v) - 2 £ £ log S^(rf «) , (8.81) 

3=1 k=l j=l 

N m 2 

• y— particles : V± = log S± v (v,j, v) — log ^(rj*" 1 " 1 ' — v) , 

j=l fc=l j=l 

N m 2 



M|fw-strings : = log Sl* v M (uj, v) - log S kM [r] 



(fc+i)- 

3=1 k=l j=l 

• M\w— strings : Vm\ w = . 

Exact Bethe equations: g& < g < g^ 

The exact Bethe equations take the following form 



N 

m(2n k + 1) = log Y u (u k ) = iLp k -J2 lo S S li(l) u k) + log (1 + *q) * 

-2£E l0 S ^(r}^- «*) + 2 log 2 - 2 £ log ( Uj - u k - |) 

fe=i j=i i=i y i fc 

+ 2 log Z M , OT * + 2 log (l - i-) * l^ 1 * + 2 log (l - -L) * i^ 1 * . (8.82) 

The consistency with the BY equations again implies the fulfillment of identities of 
the form ( |S.74 ) that we have checked numerically for the Konishi-like states. 



8.7 Reality of Y- functions 

In this appendix we show that the reality condition for Y-functions is compatible 
with the canonical TBA equations for Konishi-like states. We consider only the 
weak coupling region. The generalization to other regions and general s 1(2) states is 
straightforward. 

To start with we introduce the principle value prescription for the integrals in- 
volving log/(w) where f(u) is real for real u, and has first-order zeros (or poles) in 
the interval [a, b] at u k , k = 1, . . . , N 

log = limy] [ du\og\f(u)\K(u,v), (8.83) 

where u = a, u^ + i = b. In the cases of interest a = — oo , b = oo or a = — 2 , 6 = 2. 
Assuming for definiteness that f(u) has N zeros, and /(oo) > 0, one can write 

N 

f(u) = f(u)H(u-u k ), (8.84) 

k=l 
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where f(u) > for any uGR. 

Then the convolution terms of the form log / * K with the integration contour 
running a little bit above the real line can be written as follows 

logf*K= f du log f(u + iO)K(u + iO,v) (8.85) 

J a 

nb N rb 

= / du log / (u)K(u, v) + / du \og(u — Uk + iO)K(u + iO, v ) 

J a k=1 J a 



v) = log / *„.„ K + - log 



S(a, v) 

k=i fc=i v ' y 

where logS^w,-^ = 27ri J" du'K(u',v), and it can differ from the S-matrix defining 
the kernel K by a function of t> . It is convenient to choose the normalization S(a, v) = 
1, and most of our S-matrices satisfy this condition. 

Let us now use (|8.85| ) to show the reality of Y-functions. We start with ( |B.69| ) 



that we write as follows 
Y N 

log y(v) = - log S Uy (u v v) + log(l + Yq) * K Qy , (8.86) 
j=i 



log F + Y» = ]T log S 1 (u j -v)- log (1 + Y Q )*K Q + 2 log - - * K 
where 



| ,_ w 

3=1 Ym\w 



xJu — L Q) — x(v) x s (u + 1 Q) 7-t 

£? (o, „\ — 9 9 x ^ fS 87^1 

W ' ] " + • (8 - 87) 

and we used that 

S QV (U V) n n U-V-iQ 

S Q yM= ~ Qy ' , Sgfo - v) = S Q v(u,v)S?(u,v) = « , (8.88) 

(W,f j ) M V ~T 

and that Sq(u — v) analytically continued to the string-mirror region is equal to itself. 
Taking into account that 

2 log _Jg^ * K M = 2 log 7^ K M - ^ log S-xK - v) , (8.89) 

*M"|tu fe=l 

we see that the equations for Y± can be written as 
Y N 

log y_(v) = -Y j log ") + log(! + Yq) * K Qv , ( 8 - 90 ) 

i=i 

log K+Y_» = - log (1 + Y Q )*K Q + 2 log - | * p . v . K M . (8.91) 
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Assuming now that Yq, Y m \ vw and Y M \ W are real, the reality of Y± just follows from 
the reality and positivity of SQ tV (u,v) that can be easily checked by using ( |3.87|) . 
Next we consider Y M \ VW . By using the p. v. prescription we write ( |8.70| ) as follows 



j=l \ 3 f 



*p.». K M , M + log ^ -k K M , (8.92) 

l — 



where 



u — -(1 + M) u — -(M — 1) 
= ^ i C i . (8.94) 

Thus, we get 

t;) 2 - i) - x(« - JM) + 5) " ^^M) 



5iM (u _ v) x ^ u + i) _ ^ + i M) _ i) _ 



(8.95) 



g) 1 g) x(v-i-M) 



X 



x> - j) - a;(t; + jM) + p " ^T|m) 



Xs (u + ^)-x(v-^M) x a {u-$- 



which is obviously real and positive, and at M — 1 it has a double zero at v = u as 
it should be. 

Finally we consider lQ-functions. We write ( |S.68| ) as follows 



V 



log Y Q (v) = -LS Q -J2 (Uj, v)Sl2 x ( Uj , v) + log (1 + Y Q ,) * K%% 



' ^ "M'Q i_t nrr f-t _ ^VrsQ a™ ft ^ \^isvQ 



-lOg ( 1 + j^.^+bg (1-^7- J^+lOg ^1 - — )** 



Here 



(8.96) 



x(m - |) - s(v + |g) s(n + |) - zfo + |g) zfo - |g) 
^ *(« - J) - - fg) x(u + j) - *(« - ig) *(« + jg) 

- 1) - + |Q) 1 ~ *.(tH-iM«+iQ) 

— = M u > u )> ( 8 - 97 ) 



x s (u - i) - x(u - iQ) 1 



1 



9 y V i s (u+i)i(ii-ig) 
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where }iq(u,v) is the function that appears in the crossing relations, see [66]. Note 
that S^2x IS unitary: (S^ x ) = V^w = ^vwx/^q- 

To prove the reality of S^^(u,v)Sl^ x (u,v) we use the crossing relation. 11 To 
this end it is convenient to go to the 2-torus. Then, the crossing relation for Sj^s 
can be written in the form |69 ] 



qlQ 
D s[(2) 



( z i ~ y ^2)^ 2) ( 



UJ 2 

Y 



Z l + —, Z 2) 



1 



hq{u,v) 



2 : 



(8.98) 



and we have 



Sll(2)( Z *li Z 1, 



nlQ / W 2 



qlQ 
D S [(2) 



2 ■ 



Z 2 



(8.99) 



where we have chosen Zk to be real, and used the generalized unitarity condition for 
the mirror model. Taking into account ( |S.97| ), one concludes that S^SS" 1 ^. is real. 
It is also possible to show that the product is positive by representing S^S Sl® x as 
a product ss*, where s ~ a. 



8.8 Prom canonical to simplified TBA equations 

In this appendix we discuss the derivation of the simplified TBA equations from the 
canonical ones for Konishi-like states in the weak coupling region. Since it basically 



repeats the one done in |22|, [23[ we will be sketchy. 



Simplifying the TBA equations for v w-strings 



To derive the simplified equation (|4.3|) for fw-strings from the canonical one, one 
applies the inverse kernel (K + 1)m N to fl3.70| ) and uses the following identity 



logS^ *^ + 1)^ = 6 Ml (\ogS(u-v + --iO)+logS Uv *s > ) . (8.100) 
Note also that one should understand log 5' 1 * Q as 

log Slf{u, v) = \og{v -u)+ log S ^ ^ ' ^ . (8.101) 

v — 11 

Then, even though the last formula is valid up to a multiple of iir, it agrees exactly 
with Mathematica choice of branches. 



11 It follows closely the proof of the unitarity of the mirror S- matrix 20 



- 62 - 



Simplifying the TBA equations for Q-particles 

We want to apply (K + I^at to the equation ( |S.68|) . This requires computing 

log * (K + l)j} Q = - log S 1M ★ (K + 1)-J Q - 2 log E UM *(K+ l)j} Q 

= -5 Q2 log S(u -v)- 25 Q1 log Z u *s . (8. 102) 

As was shown in |]66|] , the improved dressing factor is a holomorphic function of the 
first argument in the string region, and the second one in the intersection of the 
region {\y%\ < 1, \y^\ > 1} with the mirror region liay^ < 0, which includes the real 
momentum line of the mirror theory. This immediately implies that 



log I^* (K + 1)^ = for Q ^ 1. 



(8.103) 



To find logS lit M*(A' + l) M 1 1 , and, therefore, Ex» we start with the kernel Kf M (u,v) 



and use the following identity derived in p3 



(8.104) 



where the kernel Kq(u, v) does not vanish for |i>| < 2, and is given by Q8.19p . Setting 
Q' = 1 = Q in ( gTOj ) and analytically continuing the first argument to the string 
region, one gets 



(8.105) 



where the analytically continued kernel Kf is given by flP3P . To obtain eq.(|8~2l) 
from QB.19Q one uses the following formula 



I Q (u+ - ,v) 
9 



h(u,v) - 



2m 



u — V 



— dt- 



u-t- 1 - 

a 



K(t,v) 



h(u,v) -K ss {u - - ,v) 



(8.106) 



Finally, integrating ( p.23 ) over the first argument, one gets (|8.24|) . 

The remaining part of the derivation of the simplified equations for Q-particles 
repeats El. 
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